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^ ■ Abstract 

In this series of articles we describe a systematic approach for studying QCD vacuum 
structure using the methods of lattice gauge theory. Our framework incorporates four 

£N) ■ major components, (i) The recently established existence of space-time order at all 

scales (fundamental structure) observed directly in typical configurations of regularized 
path-integral ensembles, (ii) The notion of scale-dependent vacuum structure (effec- 
tive structure) providing the means for representing and quantifying the influence of 
fluctuations at various scales on physical observables (phenomena), (iii) The unified 
description of gauge and fermionic aspects of the theory which facilitates a high level 
of space-time order in the path-integral ensembles, (iv) The strict "Bottom-Up" ap- 

■M ■ proach wherein the process of identifying the vacuum structure proceeds inductively 

i— h , using the information from valid lattice QCD ensembles as the only input. In this 

work we first elaborate on the meaning of the notion of space-time order in a given 
configuration which is conceptually at the heart of the path-integral approach to vac- 
uum structure. It is argued that the algorithmic complexity of binary strings associated 
with coarse-grained descriptions of the configuration provides a relevant quantitative 

^ ■ measure. The corresponding ensemble averages define the ranking of different lattice 

theories at given cutoff by the degree of space-time order generated via their dynamics. 
We then introduce the set of local transformations of a configuration, chiral orderings, 
in which the transformed gauge connection represents an effective matrix phase ac- 
quired by chiral fermion when hopping over a given link. It is proposed that chiral 
orderings facilitate the evolution in the set of actions which increases the degree of 
space-time order while preserving the physical content of the theory, and should thus 
be used in the search for the fundamental QCD vacuum structure. The relation to 
renormalization group ideas is discussed, and the first step in general formulation of 
effective lattice QCD realizing the notion of scale-dependent vacuum structure is given. 
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1 Introduction 



Understanding the vacuum structure of non-abelian gauge theories is the subject of active 
theoretical research. While no definitive answers have been obtained over a number of 
years, the subject continues to captivate the interest of many researchers since the stakes 
are high. Indeed, the potential fruits of such effort include, among other things, identifying 
the mechanism leading to confinement and spontaneous chiral symmetry breaking (SChSB) 
- phenomena believed to be central for the dynamics of strong interactions. 

While the vacuum of quantum field theory is a state in the appropriate Hilbert space, this 
is not necessarily the representation of choice in the search for vacuum structure. Indeed, the 
discussion of related issues in QCD is frequently carried out in the Euclidean path-integral 
formalism. Here one does not work with states in the Hilbert space, but rather with the 
statistical ensemble of gauge configurations, i.e. the set of pairs £Q CD = { ( A, P[A] T>A ) }, 
where A = {A^(x)} collectively denotes an arbitrary space-time configuration of gauge 
potentials, and P[A]VA oc exp(-S QCD [A])VA is its weight ("probability") in the path 
integral. Since the vacuum expectation value of any operator is replaced by the corresponding 
ensemble average, the information stored in the quantum vacuum state is equivalent to the 
information encoded in the ensemble. In this sense, one can view the ensemble as a particular 
representation of the vacuum. 

One advantage of working with classical ensembles is that the notion of "vacuum struc- 
ture" acquires a natural intuitive meaning here. Indeed, the attempts to understand QCD 
vacuum in this language aim at finding a related ensemble £ STR (ensemble with "structure"), 
which only involves a certain definite subset of configurations that share a specific kind of 
space-time order. The notion of "vacuum structure" is then loosely associated with the 
common space-time structure exhibited by these special configurations. Ideally, one should 
be able to describe £ STR analytically in terms of collective degrees of freedom (objects in the 
gauge field) that encode this order. The (usually vaguely stated) criterion for identifying the 
ensemble £ STR is that the "QCD physics is reproduced" while, at the same time, this physics 
should be naturally understood in terms of the associated order and/or the corresponding 
collective degrees of freedom. There are two highly interrelated problems with the above 
program which currently appear to prevent it from becoming a developed and systematic 
area of research. 

PI: The lack of well-posed goals consistent with the fact that £® CD describes quantum field 
theory. 

P2: The lack of guiding principles and the systematic framework for constructing £ STR once 
such goals are given. 

The aim of this series of papers is to discuss a set of ideas that have a bearing on the above 
issues. Most of these ideas are original, but some represent a significant expansion of the 
line of thought that already appeared in our previous work. In the remaining parts of this 
section, we will discuss in more detail why we view PI and P2 as pressing problems, which 
will allow us to naturally set the stage for our approach. 
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1.1 Lattice QCD 



One of the important reasons for using the Euclidean path integral formalism in the context 
of vacuum structure is that gauge theory can be conveniently regularized in this framework. 
Indeed, the considerations of the previous section are only meaningful when viewed as a 
result of the limiting procedure involving regularized versions of the theory. Gauge invariant 
regularization on the Euclidean space-time lattice [1] provides an appropriate framework for 
both the non-perturbative definition of gauge theory and the study of the vacuum structure 
along the lines described above. 

The power of lattice QCD in this context resides with the fact that the vacuum of 
a finite lattice system, namely the ensemble S L( ^ CD , is at our disposal in the following 
sense. The ensemble S LQCD — {(U, P{U)dU ) } can be fully represented by an infinite 
probabilistic chain of configurations {. . . ,U l ^ 1 ,U\U l+1 , . . .}, such that the probability of 
encountering the configuration U = { U x>ft G SU(3) } is proportional to its weight P(U)dU oc 
exp(— S L ® CD (U)) dU in the lattice path integral (a finite-dimensional integral). At the same 
time, finite sections of such chains are generated explicitly in Monte Carlo simulations of 
lattice gauge theory on digital computers [2, 3]. In other words, a numerical simulation gives 
us direct access to the vacuum in the form of configurations dominating the evaluation of 
the regularized QCD path integral which ultimately defines the theory. 

Given the above, one naturally expects lattice QCD (LQCD) to be very relevant in the 
quest for understanding the QCD vacuum structure. One of the points that we will emphasize 
throughout this series of papers is that the role of LQCD should in fact be elevated from 
auxiliary and secondary to primary and decisive. 

1.2 The "Top-Down" Approach 

For a number of years, the research in this area proceeded mainly along the line which we 
refer to as a "Top-Down" approach. The characteristic feature of Top-Down is that the 
starting point of an investigation is an immediate proposal for relevant collective degrees of 
freedom (e.g. an instantons, magnetic monopoles, center vortices). Such a proposal is made 
in the continuum and is usually supplemented by a "picture" of how confinement or SChSB 
could be understood if the collective degrees of freedom in question would indeed turn out 
to be the most relevant. 

The role of LQCD is rather auxiliary in this approach. The primary aspect amounts 
to calculations supplying details needed for an eventual analytic representation of £ STR in 
terms of the collective degrees of freedom chosen, i.e. determining the densities, sizes and 
other geometric characteristics of the objects relevant for the description of S STR . Such an 
inquiry typically proceeds via certain procedure that assigns to every configuration U of 
the numerically-generated ensemble £ L ® CD = {..., [7 1-1 , U l , U' l+1 , . . .} an associated config- 
uration Uc = Fc{U) in which the corresponding degrees of freedom C are "visible". This 
produces a numerical ensemble with "structure" £ LSTR = {. . . Uq~ , Uq, Uq' 1 , . . .} on which 
the needed characteristics of S can be measured. Another use of £ LSTR could be to 
quantify (by measuring and comparing variety of relevant observables in £ L Q CD and S LSTR ) 
how "close" £ STR is to £ QCD , but this is usually not done beyond the string tension. 

There are numerous problematic issues associated with the Top-Down approach some of 
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which we list below. 

(1) The proposal for the collective degrees of freedom is always ad hoc in the sense that it 
is motivated by certain interesting properties (mainly of topological nature) but not taking 
into account at all the global structure of quantum QCD ensemble. In that sense it amounts 
to little more than guessing (see P2). Nevertheless, it is certainly possible (but maybe not 
likely) to get to the root of a problem via a well-motivated guess. 

(2) The procedure for selecting U — > Uc = F C (U) is not unique in existing examples, and 
is sometimes gauge dependent. Moreover, the aspect of non-uniqueness of Fc appears to be 
generic. 1 The finality of the results and arguments in favor of a given picture arrived at in 
a Top-Down manner thus appears to always remain in question (see P2). 

(3) It is not clear how to interpret the meaning and the scope of S STR arrived at in a Top- 
Down manner in the context of quantum field theory defined by £ QCD . Indeed, at the face 
value S STR = £ STR (£ QCD , C, F c ) is a function of £ QCD , the collective degrees of freedom C, 
and the procedure U — > Fq{U) which makes C explicit in the configuration. It is hard to find 
a general notion in field theory which would give well-defined meaning to such generic 8 . 
To illustrate this further, it is customary in the Top-Down approach to interpret £ STR 
in terms of a phenomenon (physics) one wants to understand. For example, the ongoing 
attempts frequently start as an effort to understand confinement, i.e. to construct £ STR 
in terms of correctly chosen C such that string tension (er) would be reproduced and its 
origin understood. There are similar attempts to understand chiral symmetry breaking and 
reproducing chiral condensate. However, this approach carries ambiguities that are difficult 
to resolve conceptually. For example, having constructed £^f for confinement, do we expect 
that it will also naturally explain chiral symmetry breaking and reproduce the condensate 
(and vice versa)? Such a coincidence might not be very likely since the construction of 
^conf 1S almost certainly not unique. Indeed, it appears more reasonable to expect (and 
different coexisting pictures of confinement support this expectation) that there are many 
ways to construct S^cmf reproducing just a. Similarly, there are probably many ways to 
construct S^ R which reproduce just the chiral condensate, thus supplying multitude of 
different "pictures" for SChSB in QCD. Moreover, the argument doesn't stop there. One 
can legitimately ask how vacuum encodes (and explains) the mass of the lightest scalar 
glueball or, for that matter, the mass and properties of any of the other bound states. The 
point of these remarks is that organizing the search for vacuum structure "by phenomenon" 
is probably not the appropriate way to pose the problem (see PI). This simply reflects the 
fact that field theory is not naturally organized in terms of "phenomena" . 



1.3 New Input from Lattice QCD 

There were recent developments in this area that will serve as an important part of the 
conceptual input for our framework that we hope can overcome the issues related to PI 
and P2. To fix the notation and the language, let us denote by £^ QCD the ensemble of 

1 There is no satisfactory mathematical meaning to the task of constructing the map Fq for any reasonable 
class of C. In other words, the attempt to define the "C-content" Uc of configuration U appears to be ill- 
founded. 
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configurations = <f>(U) (local composite field) associated with given £ L Q CD . 2 In terms of 
a probabilistic chain we have 

EE {...^-\<P\^ + \...}^S^ QCD (1) 

The probability distribution (and the associated action) governing the ensemble £^ QCD [ s 
related to the fundamental LQCD distribution P(U) via 

= e~ s ^/Z = J dUP(U)6(<f>-$(U)) (2) 

We now wish to highlight two aspects put forward in recent works [4, 5, 6, 7] dedicated to 
study of a particular composite field, namely the topological charge density [8, 9] associated 
with Ginsparg- Wilson operator [10], i.e. q x = q x (U) = — tvy 5 (1 — D X;X (U)/2), if D satisfies 
the canonical Ginsparg-Wilson relation. 

II: The configurations q appearing in probabilistic chains representing £^ QCD exhibit an 
observable space-time order (structure) [5, 6]. 3 This order disappears after random- 
ization of space-time coordinates of the underlying gauge field [7] , and is thus a man- 
ifestation of space-time order in £ L Q CD . 4 

The conceptual content of the above result is non-trivial. Indeed, this is the first time that 
a definite space-time order has been observed in the physically relevant composite field eval- 
uated directly on unmodified configurations that appear in probabilistic chains representing 
gLQCD rjj-^g £ a( ^ gj veg substance to the (previously quite unfounded) expectation that QCD 
path integral is dominated by certain definite class of configurations exhibiting a common 
structure. We emphasize that this didn't necessarily have to be the case. While the ensemble 
averages build up coherence over distances of order 1 fm in QCD (as measured via properly 
defined physical correlators), it was not a priori obvious that this coherence would manifest 
itself at the configuration level via explicitly detectable space-time order. Nevertheless, it 
now appears that this is indeed so. 

12: It was proposed [4, 6] that the patterns of structure in topological charge fluctuations 
described by q associated with Ginsparg- Wilson fermions can be meaningfully studied 
as a function of the fermionic response scale A F . The underlying proposition is that 
all aspects of QCD vacuum structure should be viewed in a scale- dependent manner. 

2 Note that the symbol for configuration is distinguished from the symbol for the local map 4> relating 
it to the configuration of gauge field U. 

3 The nature of the ordered structure present in these configurations is quite intriguing (see [5, 6]) but we 
will not discuss it here. 

4 Strictly speaking, the corresponding studies have so far been performed in pure-glue LQCD defined by 
Wilson and Iwasaki gauge actions, and using the Wilson-based overlap operator [12] to define the topological 
field. However, it is expected that the existence of the structure extends to generic LQCD and to generic 
Ginsparg-Wilson operator. 
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To be more explicit, the scale dependent ensemble S^p D for q was defined as 

E LQCn = { r i q i q > + i } _^ {... r ^, q ^ F , q ^ F ,...} = E^ D (3) 

Here for arbitrary configuration q, the effective topological field q AF is defined via [4, 6] 

<£ F = - E (l-^)(^) f 75^ (4) 

A<aA F Z 

where ?/> A is the eigenmode of D with eigenvalue A and a is the lattice spacing. This formalism 
allows to study the patterns in topological charge fluctuations as a function of scale. 



1.4 Resolutions: "Bottom-Up" and the Fundamental Structure 

Given II and the host of problems with the Top-Down methodology described above, it 
appears highly desirable to dramatically change the point of view when approaching the 
problem of QCD vacuum structure in the path integral formalism. Since now there exists 
evidence for a directly observable order in typical representatives of £ ^ , it might no 
longer be necessary to rely on guessing associated with Top-Down strategy. One can instead 
take the attitude that a viable candidate for relevant collective degrees of freedom C has 
to have the space-time attributes consistent with the structure observed in typical LQCD 
configurations. Thus, instead of making fixed C a starting point of the investigation, the 
approach built here views determining C as one of the final goals in a systematic study. We 
refer to this as a "Bottom-Up" approach to QCD vacuum structure. There are two crucial 
aspects that we wish to emphasize and to associate with this approach. 

Al: The sole input used in the Bottom-Up search for £ STR comes from various valid ensem- 
bles £ L Q CD defining £^ CD in the continuum limit. At the technical level, this can also 
enter via associated ensembles £^ QCD of various local composite fields (see Eq. (1)). 

Indeed, the intention in the Bottom-Up approach is that its only source of information are 
typical configurations from £ L Q CD (for various lattice formulations), and that it starts from 
a white page, i.e. using no assumptions at all about the nature of collective degrees of 
freedom C in terms of which £ STR will ultimately be described. Consequently, this research 
will naturally proceed via intermediate phases involving the creation of suitable concepts 
and characteristics whose behavior in £ L Q CD will provide hints about the nature of £ STR . 
In particular, one envisions the following. 

(i) Geometry. The Bottom-Up approach to QCD vacuum structure should start by ex- 
ploring the geometric patterns occurring in the fields typical of £ L Q CD [5 ; H]. The simplest 
beginning is to study the structure in the basic scalar and pseudoscalar gauge invariant 
composite fields, i.e. action density and topological charge density. More complex studies of 
other composite fields, and perhaps also of fundamental gauge fields, can follow. These inves- 
tigations should lead to precise geometric statements about £ L ® CD ; which can be supported 
or disproved by numerical simulation. 
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(ii) Picture. The information from geometric studies will serve as a sole guide and input for 
an emerging picture of the relevant collective degrees of freedom. The underlying expectation 
is that the requirement that such picture be consistent with accumulated knowledge about 
geometric behavior of fields will significantly restrict the available possibilities. 

(iii) Physics. The viability of any picture consistent with geometry will ultimately have 
to be assessed on physical grounds. In other words, some relevant physical implications (or 
predictions specific for that picture) have to be compared with physics extracted directly 
from £ L Q CD , It is worth emphasizing that stages (ii) and (iii) are naturally intertwined 
in that physics considerations immediately come into focus when considering any particular 
picture of a collective variable. 

(iv) Definition of £ STR . When the qualitative and quantitative knowledge about the 
relevant collective variables in QCD reaches a sufficient level of completeness via the above 
systematic process, it might be possible to define the degrees of freedom C analytically. A 
precise definition of £ STR in terms of C then might become possible and the understanding 
of QCD vacuum could be achieved. 

A2: The vacuum structure obtained from £ L ® CD in a Bottom-Up manner is fundamental [5] 
in the sense that it involves fluctuations of elementary fields at all scales upon taking 
the continuum limit, and should thus be relevant for all aspects of QCD physics. 

Indeed, the fact that the structure (embodied in collective degrees of freedom C) is extracted 
from behavior in full unmodified ensembles £ L( ^ CD implies that, in principle, no aspect 
of physics should be completely lost in the process of identifying the structure and the 
subsequent construction of £ STR . We emphasize that this fact represents a non-trivial shift 
in the focus and the scope of the QCD vacuum research. Indeed, the line of thinking adopted 
in Refs. [5, 6] together with the results obtained indicate that there might exist a structure 
that encodes the origin of chiral symmetry breaking, confinement, the mass of the lightest 
pseudoscalar glueball 5 , the mass of the 10-th excited scalar glueball and so on - all in one 
fundamental vacuum structure. 

At the current stage, there are several missing pieces (of both conceptual and practical 
nature) that need to be resolved before the Bottom-Up approach to fundamental QCD 
vacuum structure can evolve into a developed area of research, (a) So far the practical 
framework only exists for a single (albeit highly relevant) composite field, namely topological 
charge density. While this might eventually be sufficient to uncover the nature of C, there 
should be a general framework that would allow to study an arbitrary composite field and 
the gauge field itself, (b) There is no accepted conceptual notion for degree of "structure" 
(order) in a given individual configuration of the gauge field (or a composite field for that 
matter). While the intuitive notion might be sufficient to proceed, the systematic framework 
should contain at least a concept that defines what we are attempting to do. (c) Once such 
notion is in place, there should exist a mechanism that guides the process of identifying C 
via maximizing the degree of order in £ L Q CD ensembles. This might require concepts that 

5 The masses of bound states are not always thought of as associated with vacuum structure but in very 
real sense they are. After all, the correlation functions encoding these masses are vacuum expectation values 
of appropriate operators. In fact, all the observables arise as averages over £Q CD and thus one can attempt 
to "understand" their origin in the vacuum structure. 
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do not belong to the existing framework of quantum field theory, but should nevertheless 
be meaningfully integrated into it and, hopefully, enrich it. Addressing these issues is the 
primary goal for the first two papers of this series. 

1.5 Resolutions: Scale-Dependent Vacuum (Effective Structure) 

We have emphasized in Sec. 1.2 that associating £ STR with a particular phenomenon is 
rather problematic and unnatural in the general framework of quantum field theory. On the 
other hand, it has proved very useful to organize field theory calculations (and thinking) in 
terms of a scale. Given this and 12, it is thus reasonable to expect that this will also be 
a productive way to deal with the QCD vacuum. In other words, we suggest the "change 
of basis" in which the questions about QCD vacuum are asked. In particular, instead of 
seeking ^phenomenon f° r various different phenomena, we suggest to seek S^ TR , i.e. the vacuum 
structure relevant for physics at arbitrary typical scale A. Another way of saying this is that 
the process involving momenta up to typical value A will "excite" or "see" the vacuum 
structure involving the corresponding momenta. Thus if we wish to understand the role of 
vacuum in such phenomenon, it is £'l TR that is relevant for such understanding. 

Similarly to the case of fundamental structure, we wish to emphasize two aspects asso- 
ciated with the Bottom-Up program applied to uncovering the scale-dependent structure in 
QCD vacuum. 

A3: The sole input used in the Bottom-Up search for £l TR comes from various ensembles 
gLQCD ^ ase( j^ Qn gLQCD _ Alternatively, this can also enter via associated ensembles 
£^ CD of composite fields (p (see e.g. Eq. (3)). 

A4: The structure obtained from £ A in a Bottom-Up manner is an effective structure 
[4, 5] exhibited by elementary fields dominated by fluctuations at the momentum scale 
A. The scale-dependence introduced in this manner allows to associate a typical scale 
(and the corresponding effective structure) with various phenomena. 

We emphasize that the notion of £f TR is essential for the interpretation of the fundamental 
structure represented in S STR . Indeed, the existence of such an all-encompassing structure 
raises a conceptual issue namely how does one "decode" the structure relevant for a particular 
physical phenomenon. The concept of S^ TR resolves this and represents an "unfolding" of 
the fundamental structure into an effective structure relevant for physics at scale A. Note 
that the sought for collective degrees of freedom C characterizing the fundamental structure 
can acquire a scale-dependence C = C(A) in this process. 

While the discussion above summarizes our scale-dependent approach, there are clearly 
some crucial gaps to be filled before this can be viewed as a general theoretical/practical 
framework. Indeed, so far such framework (partially) exists only for the study of topological 
charge density. How does one generalize it to other operators and to the ensemble level? In 
other words, how are E L K QCD and S^ Q A CD defined? Can they be defined for arbitrary S LQCD 
and arbitrary 0? Also, even in case of topological charge density one needs to address 
the question of how is the fermion-response scale A F (see Eq. (3)) related to the usual 
momentum scale A. We will start elaborating on some of these issues in the present work, 
and the problem will be addressed fully in the third paper of this series. 
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1.6 Final Introductory Remarks 



We wish to make three sets of comments before we start. 

(i) The ideas presented in this series of papers (and the manuscripts themselves) have been 
developed over the period of several years, and certain aspects of it were discussed by the 
author in less unified manner at various meetings. In this extended period (and before that), 
there appeared works that, at least to some extent and at least in some regards, align with 
the points of view expressed in this work. Being such, we would like to mention some of 
them. 

The idea of using low-lying fermionic modes as an unbiased probe of vacuum structure 
can be viewed as a precursor to Bottom-Up. The beginnings of the "mode" approach can be 
traced to the set of early works [13], where they were used to assess the global properties of 
the underlying configuration. The proposal for using fermionic modes for systematic study 
of local behavior in the configuration has been made in Ref . [14] . The related early papers in 
this category include [15, 16, 17, 18], and later [19]. 6 Recently, there also appeared related 
studies that use the eigenmodes of the covariant lattice Laplacian for similar purposes [20, 21]. 

As emphasized in [4], the problem with studying individual eigenmodes is that they 
do not provide a direct interpretation of vacuum properties in terms of gauge fields or 
gauge invariant composite fields. As such, they have to rely on the model assumptions 
(although these might be testable a posteriori). Focusing on studying the properties of 
gauge invariant composite fields expressible in terms of fermion eigenmodes is the natural 
way to proceed [4, 5, 6] with Bottom-Up. This approach has recently been taken up also 
in Refs. [23], which also confirmed the existence of the fundamental structure of Ref. [5]. 
Similar fundamental structure was observed in case of CP^ 1 models in Ref. [22]. Using a 
certain non-standard non-local topological density operator, vacuum properties have been 
very recently investigated in a largely Bottom Up manner in Ref. [24]. The recent work 
on using the Laplacian modes for filtering the gauge fields [21] is an interesting step in the 
Bottom-Up direction as well. 

(ii) Throughout this series of papers we will use a rather formal mathematical language 
when stating certain important conclusions. 7 We do this even though in most cases these 
statements are not rigorously proved at the time (referring to them as conjectures instead of 
theorems), and even though it is very likely that in majority of cases they might never be 
proved. The reason for this is not to do mathematics, but rather to benefit from borrowing its 
language. Indeed,the rationale for such attitude is the belief that stating intended conclusions 
precisely helps the process of obtaining meaningful results in case of investigating the QCD 
vacuum structure via path integral. The following points offer some clarification for our 
point of view. 

• If a precise statement is eventually invalidated (numerically or otherwise), we are still left 
with a non-trivial positive piece of information, namely that a particular definite possibility 

6 Note that even though the "tool" is in principle unbiased, it can be used in a biased way. Here we will 
not attempt to sort out which claims contained in all these works actually correspond to reality or can be 
interpreted as unbiased statements. 

7 We have used this approach in some of our previous works and hope to continue this practice to the 
largest extent possible also in our future "Bottom-Up" investigations. 
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is excluded. 



• A precise statement that turns out to be false can frequently be shown to be approximately 
true in a well-defined exact sense. The information value of such proposition can thus 
obviously be very large. Moreover, a false statement of this kind can frequently be modified 
in such a way (after additional information is obtained) that it captures the truth better, or it 
can even evolve into an exactly valid proposition. Indeed, we are interested in a set of precise 
statements that get better and better constrained as more numerical data is accumulated. 
Such a process has a chance to converge to meaningful definitive results. 

• In a lattice theory with finite degrees of freedom (where continuum results are obtained 
via a limiting process) everything is well defined. Moreover, as emphasized in Sec. 1.4, the 
crucial part of the vacuum structure research is of geometric nature. There is just no good 
reason not to be as specific as possible in such setting. 

(Hi) Since this work has many different aspects to discuss and our goal is to present it in a 
manner that is as coherent as possible, we decided not to include any numerical data in this 
series of papers. Rather, there will be several accompanying papers where certain conjectures 
proposed here will be supported (or possibly invalidated). Our emphasis here is on building 
of the framework indeed. 

2 The Set of Lattice Actions 

As indicated by Al, our attention will be focused on the ensembles S LQCD corresponding 
to various lattice theories S LQCD defining QCD in the continuum limit. The underlying 
idea is to build a conceptual framework which will guide us in the process of selecting 
theories judiciously in a manner that will maximize our chances of identifying the relevant 
collective degrees of freedom in QCD. To pursue this, we first need to specify the set of 
acceptable lattice actions. For most part, this discussion will be carried out in the context 
of gluodynamics (i.e. pure-glue QCD). However, we stress from the beginning that the full 
power of the approach proposed here only appears in full QCD. Indeed, certain crucial 
aspects will later be discussed in QCD with dynamical quarks. With that understanding, we 
will skip the "QCD" superscripts when denoting the actions, ensembles and related objects. 
Also, since we will be mostly working with lattice definitions of a continuum theory, we will 
mostly skip the superscript L signifying the distinction. Depending on the context, we might 
distinguish the gluonic (S G ) and fermionic (S F ) parts of full QCD action (S). 

Identifying proper lattice actions usually starts from classical expressions in the contin- 
uum and, to that end, we will use conventions where the gauge coupling only enters via the 
multiplicative factor of the gluonic part. More precisely, if A^(x) e su(3) (smooth vector 
field of traceless ant i- Hermit ian matrices), then the field strength is defined via 




(5) 



while the expression for the gauge action is 
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If tp(x) is a smooth Dirac field (tp(x) G C 12 and tp = ^ 70), 8 then the part of the action 
involving fermions reads 

S F = J d 4 xs F (x) s F (x) = ^(x)(7 MJ D M + m)^(a;) (7) 

where we assumed a single flavor for simplicity, and D^ip(x) = (<9 M + A^(x))ip(x). 

The standard construction of lattice gauge theory proceeds via introducing the hypercubic 
lattice of points n G Z 4 , and associating the gauge degrees of freedom U n ^ G SU(N) with 
links (emanating from n and ending at n+fj) while the fermionic degrees of freedom ip n G C 12 
(before quantization) are associated with sites of the lattice. We then write the candidate 
action for latticized theory with massless quarks in the form 

S G = (3 Y, O n (U, g) S F = Y $nD n , m (U, g)^ m (8) 

n n,m 

where the bare gauge coupling in the conventional form (3 = 2N/g 2 has been introduced, 
and we have allowed for the possible explicit (^-dependence of O and D. The criteria for 
selecting S G and S F are the following. 

(a) Classical limit. To establish the classical limit, the lattice structure is embedded into 
M 4 via introduction of the "naive" lattice spacing a and the relation x = an. If A^(x) 
and i/j(x) are smooth fields 9 defined on M 4 , then it is required that expressions (6,7) are 
reproduced upon transcription of (8) to the lattice via 

U n ,ix <— V exp(a J ds A^an + (1 — s)ap,fj ip n <— a 3 ^ 2 ip(an) (9) 

once the naive continuum limit a — > is taken. Here V is the path ordering symbol and \x 
a unit vector in direction \i. 10 

((3) Locality. Lattice actions are required to be local in the sense that the influence of 
distant fields on the contribution from site n decays at least exponentially in lattice units. 
In other words, dO n (U, g) / dU m:V < Ae~ a ^ m ^, with constants A and a being independent of 
U and g. 11 Similarly, the kernel D n rn (U,g) should decay at least exponentially in \n — m\, 
and \dD n>m (U, g) / dUi >u \ exponentially in \n — l\ and \m — l\. While exponential locality might 
not be necessary, it is believed to be sufficient for ensuring universality. 

8 Note that at this level we treat fermionic field as the complex-valued field and the "fermionic action" as 
the functional yielding the Dirac equation via the variational principle. We will thus view the transition to 
Grassmann variables as a part of "quantization" . 

9 By "smooth" we will understand differcntiable arbitrarily many times, even though a weaker condition 
would be sufficient. Arguments here do not depend on this detail. 

10 With the standard definition of parallel transporter and the conventional expression for lattice covariant 
derivative, the link variable transports the vector from position n + /t to n. Hence the path ordering 
prescription in Eq. (9). Simpler transcriptions for the gauge field, such as U n _^ *— cxp(a^4 Ai (na)) or U ntli <— 
exp(aA fl (na + a/2/2)) can be used as well. 

n Note that partial derivatives with respect to U niV are rather formal here since the matrix elements of 
U n v are constrained. What one really means here are the partial derivatives in the canonical representation 
of Sec. 3.4.1, i.e. dO n (u,g)/du" 
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(7) Symmetries. From the point of view of universality (and convenience) it is also desir- 
able that lattice theories preserve to the largest extent the symmetries of the action in the 
continuum. We thus require (%) gauge invariance, i.e. invariance under 

Un^ — G^E^G"^ i) n — G n i> n (10) 

with G n G SU(N); (ii) the invariance under transformations of the hypercubic lattice struc- 
ture, i.e. lattice translations and lattice rotations; (Hi) the chiral invariance of the massless 
fermionic action S F , i.e. invariance under infinitesimal Liischer [25] transformations 

i) — ► + i9^ 5 (l - RD)ip i> — ► $ + jpi9(l- DR) l5 (11) 

where R is local and [R, 75] = 0, which leads to Ginsparg- Wilson fermionic kernels satisfying 
{D,7 5 } = 2DR^y 5 D; (iv) Charge conjugation, (v) 75-hermiticity of D, i.e. = 75-D75. 

(5) No Doublers. The fermionic action S F requires additional consideration. In particular, 
it has to be ensured that the corresponding free action {U n ^ = 1) describes a single species 
of massless Dirac fermion. 12 

In what follows, we will denote by S G the set of lattice gauge actions S G satisfying 
conditions (a - 7). Similarly, the set of fermionic actions satisfying (a - 5) will be denoted 
as S F . It is worth pointing out that each element in Eq. (8) actually describes a one- 
parameter family of actions labeled by bare coupling g. While the explicit (^-dependence of 
lattice operators O and D can be present, the requirements (a - 7) have to be satisfied for 
all g. For example, the correct classical limit of S G G S G requires that 

O (U,g) — -a 4 ^trF^(0)F^(0) + 0(a 5 ) (12) 

under transcription (9) independently of g. 

Finally, we wish to remark on the conventions that we will follow when denoting lattice 
entities. The distinction between lattice coordinates n,m, . . . and the physical coordinates 
in M 4 , such as x = an,y = am, ... is made explicit by the change of notation. However, 
some confusion can arise due to the fact that while everything in Eq. (8) is dimensionless, we 
sometimes wish to write lattice relations in the continuum-like manner where the appropriate 
dimensions are in place. In that regard, we will follow the trend in the literature and not 
distinguish the symbols for dimensionless and continuum- like (possibly dimensionfull) lattice 
entities. However, to help distinguish the two cases, we will place the space-time index as 
a subscript if the entity is viewed as dimensionless, while we put it in parentheses if proper 
dimension is implied. Thus, for example, we have D nm = D xy = aD(x,y) = aD(n,m), 
with the relations x = an, y = am implicitly understood. Similarly, the action S F of 
Eq. (8) can be equivalently written as J2 x , y i>xD x , y (U, g)ip y or a 4 J2 n ,m i>(n)D(n, m; U, g)ip(m) 
or a A J2x,y4>(. x )D(x,y;U, g)ip(y). Note also that the link variable (being dimensionless) can 
take all possible forms, i.e. U n ^ = U x ^ = U(n,/i) = U(x,/i). 

12 It is worth pointing out here that the conditions we stated admit a highly singular behavior in the regions 
of configuration space far from classical behavior. These might become relevant for example if one claims 
that the corresponding classical ensemble (after quantization) can be fully represented by the probabilistic 
Markov chain. We will not deal with associated subtleties here. 
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3 Kolmogorov Entropy 



In this section we start building the framework for a systematic exploration of a fundamental 
vacuum structure using lattice QCD in a Bottom-Up manner. Relying on II we will accept 
the proposition that even though the definition of QCD in the path integral formalism 
involves inherent randomness, there exists a non-trivial space-time order in configurations 
that dominate the evaluation of the associated regularized path integral. In an attempt to 
uncover the nature of this order, and thus the nature of presumably existing underlying 
collective degrees of freedom C, we wish to study lattice definitions of QCD which tend to 
maximize the order and minimize the randomness. The goal of this section is to associate 
the amount of structure (order) in typical configurations with their Kolmogorov entropy, 
which we define as the Kolmogorov complexity (algorithmic complexity), of binary strings 
representing coarse-grained configurations. 13 

The need for such conceptual step is quite clear. Indeed, while the notion of order 
(structure) is apparently central to the study of QCD vacuum via Euclidean path integrals 
(certainly so in the Bottom- Up strategy), its precise meaning is usually left unspecified. 
There are at least two reasons for this. First, it might be assumed that, intuitively, it is 
quite clear what we mean by saying that some configuration is more ordered than other. This 
is certainly true in extreme cases, e.g. when comparing the configuration with constant field 
strength to a configuration with randomly generated links, but the question becomes fuzzy 
in less extreme cases. Secondly, the notion is highly non-trivial conceptually, and touches 
upon the involved issues related to proper definition of "randomness" . We wish to emphasize 
in this regard, that our interest in these issues at this point is strictly theoretical. We will 
point out that there exists a well-defined measure of space-time order in any configuration 
of lattice gauge field (and various composite fields), which allows to quantify the degree of 
space-time order for any lattice theory. In other words, we will discuss in what sense it is 
meaningful to say that theory Si exhibits more space time order than theory S2, but will 
not attempt to make these concepts into a computational scheme. 

In statistical mechanics, the discussion of "order" usually goes hand in hand with the 
concept of entropy. Consider, for example, two classical statistical systems consisting of 
identical degrees of freedom, but involving different interactions among them. We are inclined 
to say that the system with lower entropy is more ordered than the system with higher 
entropy. In the same spirit, consider the pure-glue lattice gauge theory S(g) G S G . Upon 
quantization, it becomes a classical statistical system with the associated canonical ensemble 
£(g) : where the action plays the role of "energy". We could thus use the corresponding 
Gibbsian entropy of £(g) as a quantitative measure of "order" associated with that theory 
(at any fixed g and fixed lattice volume). The problem with such definition is that it obscures 
(and avoids) the notion of the order at the configuration level - the concept we are primarily 
interested in when we try to identify the collective degrees of freedom of QCD in a Bottom- 
Up manner. In other words, Gibbsian entropy is not assigned to a given configuration 
independently of an ensemble we are studying: Gibbsian entropy is strictly an ensemble 
concept. 

13 We should point out that the term "Kolmogorov entropy" is used, with different meaning, also in the 
theory of dynamical systems (metric entropy) . However, since the measure discussed here is indeed entropy- 
like, we find it appropriate to use it as well. 



13 



Interestingly, the issues with similar ingredients arose in the information theory few 
decades ago. In particular, the Shannon's classical theory of communication [26] treats the 
objects to be communicated as outcomes of a random source with associated probability dis- 
tribution. The central notion here is the concept of Shannon entropy which only depends on 
the distribution itself, and does not explicitly refer to the "amount of information" associated 
with an individual object. The dual point of view was developed later by Solomonoff [27], 
Kolmogorov [28] and Chaitin [29], and is known as the algorithmic information theory, Kol- 
mogorov complexity, or algorithmic complexity. Here the focus is on the object itself and it 
defines the information (complexity) associated with an object with the smallest number of 
bits required to describe it via effective computation. The standard textbook exposition to 
this subject is Ref. [30], and the earliest work we are aware of on utilizing this framework 
in physics is Ref. [31]. Bellow we discuss in what sense we can carry these concepts into the 
study of QCD vacuum structure. Since the notions of information theory needed are not 
standard in this area, we will describe them in some detail before dealing with field theories. 

3.1 Shannon Entropy 

Suppose we have a finite or countable set of objects {source words) iGl, which are intended 
to be communicated to the receiver by the means of a binary string. Enumerating these 
objects in an arbitrary fixed manner, we associate each x n with its order number n (for the 
purposes of communication) which, in turn, can be mapped surjectively to the set of finite 
binary strings 23* = {0, 1}*, e.g. via canonical lexicographic assignment 

<-> e, 1 <-> 0, 2 <-> 1, 3 <-> 00, 4 <-> 01, 5 <-> 10, ... (13) 

where e is an "empty word". Given the maps, we will interchangeably denote by x the 
object, the integer or the string. 14 

For a meaningful communication, the sender and the receiver have to agree both on the 
above maps, and on the method which reproduces x from its code y. This is characterized by 
a definite decoding function D : 23* — > 23* such that D(y) = x. While D does not necessarily 
have to have an inverse (there could be several code words representing the same source 
word), in this discussion we assume that the encoding substitution E(x) = {y : D(y) = x } 
contains a unique element y. If l(z) denotes the length (in bits) of arbitrary string z, then 
the efficiency of coding a particular x = D(y) is given by l{y) = L D (x). 

When communicating a stream of source words, it is advantageous to consider only 
codes for which recognizing the end of the current code word in the stream does not require 
observation of the subsequent words. Such codes are referred to as instantaneous codes or 
prefix codes. Formally, the subset A C 23* is prefix-free, if xy ^ A (concatenation of x and 
y is not in A) for any non-empty x, y G A. The code specified by D is then called a prefix 
code if its domain is prefix-free. 

Assume that the source words are generated by a random source specified by the prob- 

14 We emphasize that especially in the case of integer versus associated binary string, we will mostly make 
no distinction (verbal or otherwise) whatsoever, and use e.g. x in the context of an integer even if it was 
previously referred to as a string. 
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ability distribution P(x) on X. 15 In other words, we are interested in the problem of 
communicating the class of messages, such that the source words are distributed according 
to P(x). The effectiveness of the communication accomplished via particular decoding func- 
tion D is then measured by the average code-word length, i.e. = ^2 x P(x) Ld(x). A 
natural problem in this setting is to ask what is the minimal possible average code-word 
length, i.e. L = min D {L D }, where the minimization goes over all one-to-one prefix codes 
D. The result is the Shannon's Noiseless Coding Theorem [26, 30]. 

Theorem 1 If H(P) = — J2 x P( x )^ n P( x ) ^ s the Shannon entropy associated with P, then 

H{P) < L < H(P) + 1 (14) 

The optimal average code-word length is thus (practically) equal to the Shannon entropy. 
The above discussion attempted to highlight two facts. (1) The concept of Shannon entropy 
is completely independent of the nature of the objects from X that source-words represent. 
In "information" terms, H(P) is the amount of information the observer (aware of P) gains 
when he witnesses an outcome of P. It is a measure of uncertainty in P. (2) We are 
led to Shannon entropy via the minimization of the "cost of description" of the ensemble 
{ (x, P(x)) }. Indeed, the message (binary string) y such that D(y) = x can be viewed as a 
description of the binary string x representing some real object from X. However, the optimal 
coding of x (its optimal description) is completely determined by P(x) and not x itself. 

3.2 Kolmogorov Complexity 

The notion of Kolmogorov complexity represents an attempt to quantify the amount of infor- 
mation associated with any particular binary string x. If such string completely characterizes 
some real object from X (i.e. it is not just a string assigned to it by ad hoc enumeration 
procedure, but faithfully reflects all of its attributes and no more), then it would provide an 
absolute amount of information stored in that object. 

It turns out that a concept can indeed be defined in such a way that it will satisfy the 
above features for asymptotically long strings. To explain the basic idea, we consider the 
same setup as in the previous section in that we consider the set of decoding functions D: 

— > N° (or equivalently 03* — > 03*), and call y such that D(y) = x a description of x under 
D. 16 We do not need to assume that D is one-to-one (there can be several descriptions of 
the same object after all) and, to have maximal generality, we initially do not require the 
domain of D to be prefix-free. In what follows, we will refer to this general set of functions 
as partial functions, where the term emphasizes that a particular function can be undefined 
for some arguments y. Even though our goal is to define the amount of information in the 
individual string x (quantified by the length of the minimal description y), it is obvious that 
we need to consider x in the context of all other strings from 03*. Indeed, in the opposite 
case we could simply choose any decoding function D such that D(l) = x, and thus associate 

15 This can clearly be thought of as a distribution on 03* with the corresponding mappings assumed, and 
with the probability of "unused" strings occuring being zero. 

16 We will follow the convention where the set of positive integers (natural numbers) is denoted as N, while 
the set of non-negative integers as N° . The associated finite subsets of consecutive integers will be denoted 
as {1,2,..., m} = N m and {0, 1, . . . , to - 1} = N^. 
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a single bit of information with any particular x, which is clearly not satisfactory. Thus, a 
natural way to proceed would appear to try to identify a function D a such that the associated 
minimal description is optimal for each x. More precisely, 

L Do (x) < L D {x) Wx,WD L D (x)=mm{l{y) : D(y) = x} (15) 

where we assume that if for given D there is no y such that D(y) = x, then Le>(x) = oo. 
However, it is easy to see that the partial function D with above properties doesn't exist. 
Indeed, assume the opposite and choose xi, x 2 such that L Do (x 1 ) < L Do (x 2 ), with the 
corresponding minimal descriptions x\, x\ (i.e. D (x\) = X\ and D {x^} = x 2 ). Define the 
function D' {y) = D {y) for all y ^ x*,x 2 , and D' {x\) = x 2 , D' (x 2 ) = x 1 . This obviously 
leads to shortening of the minimal description for x 2 and, in particular, L D i g {x 2 ) = L Do (xi) < 
Ld (x 2 ) thus producing a contradiction with the premise that D Q is optimal. 

The above argument shows that, in order to get a meaningful concept, it might be fruitful 
to relax the notion of an "optimal" function. 17 A natural proposal in this direction is to 
consider additive optimality. In particular, the partial function $ is additively optimal (uni- 
versal) for the set of partial functions if for every D there is a constant 7($, D) (independent 
of x), such that 

L$(x) < L D (x) + 7($,£>) Va;,VD (16) 

Assuming that a universal function exists, there is no reason for it to be unique. However, the 
optimal lengths under two universal descriptions clearly differ from each other at most by an 
x-independent constant, i.e. iL&^x) — L$ 2 (x)| < 7(^1, $2)- While this means that L$ will 
not define a strictly absolute amount of information in a finite string, there is universality here 
which applies for asymptotically large ones. To see that, we first point out that any L^(x) 
assigning finite value to all finite x (L^(x) obviously has to be of this kind), has an unbounded 
lower envelope L 1 g m (x) = mm z > x {L D (z)}. Indeed, let *B(/) denote the set of all binary strings 
up to length I. This is a finite set, and thus there exists x(x) = max{D(y), y e ?8(L™ n (x))}. 
By construction, for all z > x(x), we have L r g m (z) > L™ n (x), since all possible encodings of 
length L™ n (x) are simply exhausted for such z. Repeating this argument for x(x) + 1 and 
recursively shows that L 1 g m (x) (and thus L D (x)) is unbounded. For two universal functions 
$1, $ 2 w e then have 

< — — -— — >■ for l(x) — >■ 00 (17) 



L^{x) L^x) 

In this sense the additive optimality is universal for asymptotically large strings. 

The above discussion suggests that one could define the amount of information in x as 
L$(x) using arbitrary (but fixed) $ from the equivalence class of universal partial functions. 
However, it turns out that there is still no additively optimal element $ on the set of 
partial functions. Indeed, assume the opposite and consider arbitrary increasing sequence 
of universal minimal encodings x\ < x 2 < x\ . . . (i.e. = xi, ^(x^) = x 2 , ...). 18 

17 Indeed, it is easy to see that definition (15) would be meaningless even if we were interested in describing 
x in the context of finite number of other strings, let alone in the context of 53*. 

18 Note that the minimal encoding substitution Ed(x) = {y : D(y) — x, l(y) = Ljj(x)} can be made into 
a one-to-one function by fixing y = x* that appears first in the lexicographic ordering of binary strings. 
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Choose an arbitrary subsequence z* = such that log 2 x* < (log 2 x*^)/2, and define the 
new function $' via §'(y) = &{y) for y 7^ £*, and ^'(^iO = ®( x *j(i))- By construction, the 
encodings under $' for all will be shorter than encodings under $ and, in particular, 
L§i{xj(i)) = l(x*) < l{x*u\)/2 = L<s>{xj(i))/2. Here the inequality follows from the fact that 
in the canonical correspondence (13) between integers and binary strings we have l(x) = 
[log 2 (a; + 1)J, where |_-J denotes the floor function. We thus arrived at the contradiction 
since L$/ is not related to L$ via ^-independent constant. Consequently, the set of partial 
functions does not have the lowest element under additive optimality. 

While this appears to be a serious blow for the development of the concept we are 
interested in, it actually turns out to be a virtue in the end (see bellow). To expose the 
basic trick, let us consider arbitrary enumerable subset A = {Di, D 2 , D 3 , . . .} of partial 
functions. We can construct the partial function (specification method) D (y) that minorizes 
all elements of A under additive optimality in the following way. Fix certain one-to-one 
correspondence (n, y) < — > z = (n, y) (i.e. bijection N° x N° < — > N°). If string (integer) z 
represents the pair (n,y), i.e. z = (n,y), then D Q (z) = D n (y). By construction, the function 
D Q can emulate all the functions D n by properly selecting its arguments z. Let us consider 
the "description length" function associated with D Q) i.e. 

L Do (x) =mm{l(z) : D (z) = x} = mm{l((n,y)) : D n (y) = x} (18) 

and compare it to Lo k for arbitrary (but fixed) k. Denoting by x* D the canonical minimal 
encoding of x under D, there exists x-dependent i such that 

L Do {x) = l((i,x* D .)) < l((k,x* Dk )) = l(x* Dk ) + c k = L Dk (x) + c k (19) 

where c k = l({k,x* Dk )) — l(x* Dk ) does not depend on x if the pairing bijection (n, y) < — >■ z = 
(n, y) is chosen suitably. For example, with one of the standard prescriptions (n, y) = ny 
(string n concatenated with string corresponding to y), where the string n = l n 0n (meaning 
a string of n symbols 1 followed by followed by binary string corresponding to integer n), 
we have c k = 21 (k) + 1. Thus, indeed, the function D Q minorizes all D k e A. Assuming 
that things can be arranged so that D itself is an element of A, then D Q = $ is a universal 
element for A under additive optimality. 

The notion of Kolmogorov complexity is built on two basic results from the theory of 
computation, namely the fact that the set of recursive (computable) functions is enumerable, 
and that the minorizing functions constructed according to (18) are themselves recursive, 
thus representing universal elements. The fact that the set of recursive functions allows for 
a meaningful definition of minimal description length based on additive optimality lends the 
feeling of inherent "rightness" for the concept attempting to quantify the absolute amount 
of information stored in a string. Indeed, the use of a decoding function D in communi- 
cating an information on object x might be viewed as pointless unless the message can be 
effectively decoded by computing machines. At the same time, it is generally accepted that 
the notion of "computability" (or effective computation) is equivalent with the notion of 
recursive functions, i.e. functions that can be implemented on a Turing machine. 19 

19 More precisely, the partial function D{y) = x (from N° to N°) is recursive if there is a Turing machine 
T such that for inputs y for which D(y) is defined, T outputs D(y) in a finite number of steps, while it runs 
forever on inputs for which D{y) is not defined. 
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Focusing on recursive functions, the language of information theory starts to mesh with 
the language of computation theory. In the former case the string y such that D(y) = x is 
referred to as a code for x under D (or simply an argument of function D), while in the latter 
case one refers to the input string y placed on a tape of a Turing machine T implementing 
D as a program producing the output x. The enumerability of all Turing machines (and 
thus of partial recursive functions) can be easily seen from the fact that one can enumerate 
all possible finite lists of rules that define the function of a valid Turing machine. In fact, 
each list of rules can be mapped into a unique binary string. We will assume the canonical 
construction of this type, which is straightforward but tedious, and can be found e.g. in 
Ref. [30]. Searching in the lexicographic order through all strings, let n(i) be the i-th string 
encoding a possible Turing machine, which we denote as T; L . This enumeration of machines 
induces an associated enumeration of recursive functions Dj implemented by Tj. 20 The 
crucial point is that the canonical enumeration is recursive in the sense that there is a Turing 
machine which can decide whether a given string n encodes a valid machine. From this it 
can be seen rather easily that applying a construction (18) to sequence {D^i = 1,2,...} 
yields a minorizing D Q which is itself recursive. Indeed, one can construct the machine 
T u which takes the pair (i,y) <-> (i,y) as its input, and will simulate (on its own tape) 
machine Tj running program y. Given i, this machine will start generating binary strings in 
lexicographic order and check whether they represent the list of rules for a Turing machine. 
When encountering i-th such string n(i), it will run y using the set of rules specified by n(i) 
it found. By construction the machine T u computes D Q . 21 

To summarize, what we have described above are the basic logical ingredients leading 
to the fact that the set of partial recursive functions contains a universal element $ = D Q 
under additive optimality. This element has been identified as a function computed by the 
universal Turing machine T u , i.e. machine that can simulate the actions of any other Turing 
machine. Thus, the length of the shortest description under $ (namely L$(x)) provides an 
absolute measure for information stored in an asymptotically long string x. It is the number 
of bits that have to be transmitted to the receiver who can decode it in an effective manner 
(using a Turing machine). In the language of computation theory it is the length of a shortest 
program that, without any additional input, outputs x and halts when run on T u . 

We wish to make few comments at this point: 

(1) There are infinitely many universal Turing machines and infinitely many universal func- 
tions. Due to the universal nature of the concept discussed above, we can fix the universal 
machine and the universal function to be the canonical elements described above. Then the 
plain Kolmogorov complexity of string (natural number) x is defined as 

C(x) = L$(x) $ computed by reference machine T u (20) 

(2) Kolmogorov complexity formalizes the intuitive notion that strings with lots of regular 
patterns can be described succinctly (e.g. the string 11111 . . .1 ), while strings without any 

20 It should be emphasized that while T i} Tj always represent distinct machines for i ^ j, they might be 
implementing the same function. The set of recursive functions is defined as the set of all distinct elements 
in the sequence {Di, i = 1,2,. . .}. 

21 We note that the above description of T u -s function is not the proof of its existence. Such proof has to 
demonstrate the existence of the set of rules (list of quadruples of integers) defining the machine that will 
function as described. 
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appreciable regularities are not amenable to short description no matter what method is used. 
Thus, apart from information-theoretic meaning of the unavoidable amount of information 
to be transmitted, 22 it is also an absolute measure of order present in the asymptotically 
long string. In particular, when considering strings X\ and X2 of the same length, then we 
view xi as more ordered than X2 if C(x\) < C{x2). 

(3) Contrary to the discussion of Shannon entropy, we have not restricted ourselves to the 
prefix (instantaneous) codes D(y) so far. However, it turns out that the whole logical con- 
struction described above can be repeated if one considers only partial recursive functions 
with prefix-free domains. Such functions are implemented by prefix Turing machines, which 
only yield a halting calculation on the elements of such domain. Fixing a canonical enumer- 
ation of prefix machines 23 one can construct the universal prefix machine T u according to 
recipe described above and use it to define the universal prefix function \1/ and the corre- 
sponding prefix Kolmogorov complexity 

K(x) = Lq,(x) ^ computed by reference prefix machine T u (21) 

The functions C(x) and K(x) are asymptotically equal. However, K(x) has several desirable 
properties which C(x) does not. 24 Because of this (and also to continue the comparison to 
Shannon's theory) we will use the prefix version of Kolmogorov complexity in what follows. 

(4) While the concept of Kolmogorov complexity is built on the notion of effective com- 
putation, the complexity functions C(x), K(x) are themselves not effectively computable 
(recursive). Thus, the significance of the framework is mainly conceptual and theoretical. 
Nevertheless, many of its aspects have been successfully utilized in many varied areas. Also, 
the basic (and apparently deep) idea of minimal description length frequently serves as a 
guide for devising less general (and typically less elegant) notions that can be applied in 
practice in a more quantitative manner. 

(5) An elementary result from the theory of Kolmogorov complexity that will be useful for 
us is that, if f(x) is an invertible recursive function, then the complexities C(x) and C(f(x)) 
differ by at most a constant independent of x. This is a simple consequence of the fact 
that to arrange for evaluation of f(x) on an appropriate universal machine, one only needs 
to specify (input) the ordering number i of the machine Tj that computes /. Analogous 
statement obviously holds for K(x) as well. 

(6) It is sometimes instructive and useful to think about Kolmogorov complexity in terms 
of programs written in high-level universal languages (such as C ++ , LISP or Java). While 
the whole approach can be formalized and even used for practical purposes, we will not be 
concerned with this aspect here, and view it only as a tool for developing intuition. 

22 We note in this context that strings for which C(x) < l(x) (i.e. they allow shorter description than 
themselves) are referred to as compressible, while the rest of them as incompressible. 

23 One way to proceed is to use a standard procedure of transforming any given Turing machine into a 
prefix machine. This way the canonical enumeration of Turing machines directly induces the enumeration 
of prefix machines. 

24 0ne relevant difference is that, unlike plain complexity, the prefix version is subadditive, i.e. K(xy) < 
K(x) + K{y) + 0(1) (here xy means concatenation of strings x and y). 
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3.3 Kolmogorov Entropy: Z2 Gauge Theories 

The preliminaries in the previous two subsections allow us to come back to the question 
of how to quantify the level of space-time order in lattice gauge theories in an interesting 
manner. Before dealing with SU(N) theories, it is first useful to take a simpler case. For 
example, consider the set of pure gauge theories S G with N=2. To every S G e S G assign 
a Z 2 gauge theory S GZ by replacing SU (2) variables by Z 2 variables, and the path integral 
by a corresponding discrete sum. This defines the set S GZ of Z 2 gauge theories. Our goal 
is to answer (in principle) the following generic question. Take two theories S GZ (p 1 ) and 
S GZ (/3 2 ) (fls fixed) both defined on a finite discretized symmetric torus, i.e. both involving 
L 4 lattice sites. Which theory possesses more space-time order? 

Our aim is to approach this issue via concept that allows for even more detailed question. 
In particular, given two Z 2 configurations Ui, U 2 , which one exhibits more space-time order? 
We propose to answer this question using the notion of Kolmogorov complexity. Indeed, for 
fixed L, the set X of objects we are interested in is the set of all possible Z 2 configurations 
U. It is clear that every U G X can be completely described by a binary string z(U) of length 
4L 4 . Indeed, the two possible values of the field U(n,fj) associated with link (n, jj) can be 
represented by the two binary values b G 23 = {0, 1}. Moreover, if we fix some one-to-one 
function F(n, [/,) — k G {1, 2, . . . , 4L 4 }, then we can form a corresponding binary string via 

z F (U) = hb 2 . . . 6 4L4 b k = U(F-\k)) (22) 

We need to be more specific about the nature of the function F, which has to be defined for 
all L. In other words, L has to be one of the arguments of F, and also (a dummy) component 
of the image so that F is formally invertible. To do that, let us first agree that the coordinates 
n = (ni, n 2 , 713, n 4 ) on the torus of size L take the values n« G {0, 1, . . . , L — 1} = and 
remind ourselves that ji G {1, 2, 3, 4} = N4. Then we consider invertible partial functions 

F : N x (N ) 4 xNihNxN such that (N^) 4 x|h N 4L 4 under F(L,n,fi) (23) 

for arbitrary fixed L, where we skipped the L-component of the domain/range in the second 
part of the expression (see also equation (25) bellow). We then define the Kolmogorov entropy 
of configuration U relative to F via 

E*(U) = K(z F (U)) (24) 

i.e. we associate it with prefix Kolmogorov complexity of zf(U). 

The degree of order detected in z F (U) via its prefix complexity is a result of two factors. 
The first is the level of inherent space-time order present in configuration U, while the second 
is the influence of map F which, although fixed for all U, appears rather arbitrary at this 
point. However, we will argue bellow that this arbitrariness is innocuous for asymptotically 
large configurations U if we restrict ourselves to recursive functions F. Indeed, denoting the 
set of invertible recursive maps satisfying (23) as J 7 , let us consider Fi, F 2 G J 7 , and compare 
the complexities K(z Fl (U)), K(zf 2 (U)). Since F 1: F 2 are recursive, the binary strings zf^U) 
and z F2 {U) are recursively related via invertible map 

z F2 (U) = G F ^(z Fl (U)) VU where (G F ^(x)) t ee x Fi{F -i {i)) (25) 
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and we have denoted the i-th bit of string x by X{. Consequently (see comment (5) in previous 
section), the complexities K(zf 1 (U)), K(zp 2 {U)) differ at most by an [/-independent constant 
7(^1,^2). We can thus conclude that Kolmogorov entropy (24) is additively optimal with 
respect to J 7 , with arbitrary F G T being a universal element. 

It should be remarked that since additive optimality is inherent in the notion of Kol- 
mogorov complexity already, it is only natural that the Kolmogorov entropy (24) has an 
additive arbitrariness related to the choice of F. While we can choose and fix any F G J 7 , 
it is relevant here that the complexity of F can affect significantly the value of Kolmogorov 
entropy in small configurations, and possibly mask the true level of their space-time order. 
Thus, contrary to the selection of universal machine, there is a natural general preference 
for the choice of maps F, namely for those with low complexity. In fact, given T u , the 
appropriate choice F u can be constructed as follows. Let n(i) be the sequence of natural 
numbers such that T n u\ is the i-th Turing machine computing function G T . The set 
{K(n(i))} is bounded from bellow and thus has a minimum K min . We define 

F u = F k such that k = min{ % : K(n(i)) = K min } (26) 

In other words, F u is specified by lexicographically first Turing machine computing an ele- 
ment of JF, whose complexity is equal to K min . 25 Similarly to fixing universal machine T u , 
we now fix the map F to F u and thus define the Kolmogorov entropy of configuration U as 

E K {U) = K(z(U)) z(U) computed via F u (27) 

Given two configurations U\, U 2 on the same space-time lattice, we define U\ to exhibit more 
space-time order than U 2 if E K (Ui) < E K {U 2 ). 

We are now in the position to quantify the degree of space-time order associated with 
any theory S = S(f3) G S GZ . In particular, we identify this measure with the average 
Kolmogorov entropy in ensemble S associated with S, namely 

E K [S] = {E K {U))s = Y,E K {U)Ps{U) P S (U) = e~ s ^/Z (28) 

u 

The meaning of E K [S\ is quite clear. If the ensemble S is dominated by configurations 
with high level of space-time order, then its Kolmogorov entropy is low. Equivalently, if the 
individual configurations typically appearing in probabilistic chains representing S exhibit a 
high degree of space-time order, then the Kolmogorov entropy of the theory is low. 

We wish to make few important remarks at this point: 

(1) For fixed finite L, Kolmogorov entropy E K [S,L] induces an order relation on the set of 
actions S . 26 In particular, we say that S\ < S 2 if E K [S\,L\ < E K [S 2 ,L]. Since the 
elements of S GZ are local, it is expected that any two actions can be consistently compared 
in the infinite volume limit, i.e. that there exists L 12 such that E K [Si,L] < E K [S 2 ,L] for 

25 It is worth emphasizing that one defines the complexity of a general recursive function / as K(f) = 
mh\i{K{i) : Ti computes /}. The reason is that, assuming a canonical order of machines, specifying i is all 
that is needed to specify the map. 

26 More precisely, it induces a linear order on the set of equivalence classes, where Si ~ S2 if E K [S{\ = 
E K [S 2 ]. 
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all L > Li 2 . As discussed extensively above, the nature of Kolmogorov entropy is such 
that there is an inherent arbitrariness by a constant (related to the choice of universal 
machine and the choice of map F) , and the notion acquires an absolute meaning in the limit 
L — > oo. At the same time, the above ordering relation is expected to be affected very little 
by this arbitrariness even at finite L since T u and F u are fixed for all U and all theories. 
Nevertheless, if one wants to ensure strict universality in the infinite-volume limit, then 
bounded differences of entropies need to be considered as irrelevant. In particular, if there 
is an L-independent constant C 12 such that 

\E K [S 1 ,L]-E K [S 2 ,L]\ < Ci 2 VL (29) 

then Si and 62 should be viewed as possessing the same amount of space-time order. In 
other words, it is natural in the L — > 00 limit to consider the equivalence classes of theories, 
with the equivalence relation defined by (29). 

(2) In definition (27) we took into account that it is desirable (but not necessary) to minimize 
the influence of the map F on the complexity K(z(U)). This can be done alternatively 
(and in some ways more elegantly) by considering the conditional Kolmogorov complexity 
K(z(U)\F). Roughly speaking, this means that information (and necessary input into T u ) 
about F doesn't count as a contribution to the complexity of z(U). We will not discuss this 
possibility here in detail. 

(3) Since we are dealing with gauge theories, it is natural to ask how gauge freedom enters 
the issue of inherent space-time order. In particular, given gauge transformation g, the 
Kolmogorov entropies E K (U) and E K (U 9 ) can be quite different. In other words, E K (U) 
is not gauge invariant in the usual sense (at the configuration level). We can see three 
acceptable views concerning this issue, (i) Gauge invariance is not really necessary for the 
concept characterizing the space-time order in the configuration. 27 Indeed, it is simply 
a fact that gauge freedom can introduce an unphysical space-time noise. Since this noise 
is integrated over in the path integral for all theories in the same manner, the ordering 
of theories via E K [S] preserves its intended meaning even if one is interested in gauge- 
invariant quantities only, (ii) We can make E K (U) gauge invariant and define the associated 
Kolmogorov entropy for the theory via 

E K (U) = minK(z(U 9 )) E K [S] = (E K (U)} S (30) 

In other words, the gauge invariant Kolmogorov entropy of U is taken to be the minimal 
Kolmogorov complexity of z(U 9 ) over the gauge orbit {U 9 } of U. Since K(x) is bounded 
from bellow, the minimum indeed exists for arbitrary L. (iii) Instead of taking descriptional 
complexities of z(U) as a measure, we could use the complexities of strings z((f>(U)) describing 
some fixed reference gauge invariant composite field 4>(U). This possibility is less general, 
but can certainly be exploited for comparing the theories (via E^(U)) with respect to a 
given relevant composite field (i.e. topological charge density in case of QCD). 

27 Note that considerations related to Elitzur's theorem do not play a role here since E K {U) is a global 
quantity for which no sensible local definition exists. 
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3.3.1 Information Issues and the Information Gauge 

It is interesting to elaborate on the information theory viewpoint of the definitions dis- 
cussed above. Assume that the sender wishes to communicate a complete information on 
Z2 configuration U to a receiver via message decodable using effective computation. Then 
the (additively) optimal amount of information necessary to transmit is given by E K {U). 
However, if we are interested (as usually is the case) in the gauge-invariant content of the 
configuration, then we can further minimize the cost of transmission by communicating only 
E K (U) bits to the receiver. In effect, instead of using U, we are sending a description of 
a particular configuration V = V(U) from the associated gauge orbit, which possesses the 
highest level of space-time order. Such V(U) may not be unique, but we can make it unique 
e.g. by selecting a configuration which appears first in the lexicographic order of associated 
binary strings. More precisely, if V{U) is the ordered set of such configurations, namely 

V(U) = {U 9 * : E K (U 9 *) = E K {U), z{U 9l+1 )> z(U 9 *)} then V(U) = U 91 (31) 

where the ordering of binary strings is defined by the canonical correspondence (13). The 
above prescription assigns a unique configuration to a gauge orbit of arbitrary U, and thus 
can be viewed as a gauge-fixing prescription. 

Assume now that the sender intends to communicate to the receiver information about the 
gauge-invariant content of the Z2 gauge theory with action S, and he wishes to do that via 
transmitting the description of some relevant configurations. This sender would generate N 
independent configurations Uj from distribution Ps{U) and, instead of transmitting minimal 
description of Uj, he would transmit the description of Vj fixed by prescription (31) in order 
to minimize the cost. In other words, he would prefer "fixing the gauge" accordingly. For 
this reason, we refer to the gauge specified above as the information gauge. 

It is also important to emphasize at this point that the role of information theory in 
our construction is purely to assess the level of space-time order exhibited by the configura- 
tion (theory), and we are not attempting to connect it to the absolute amount of "relevant 
information" stored in it. The basic concepts characterizing the latter, though of utmost 
interest, are yet to be developed. In fact, one should readily appreciate the apparent contra- 
dictory tendencies in that regard. Indeed, the configurations with high Kolmogorov entropy 
require large number of bits to describe, and in that sense one can view them as storing a 
lot of information. However, by the very nature of our approach, we expect their relevance 
for physics to be rather limited. On the other hand, we hope that configurations with low 
Kolmogorov entropy, exhibiting a specific kind of space-time order and allowing for a short 
description, will contain large amount of the physically relevant information. 

3.3.2 Shannon Entropy and Kolmogorov Entropy 

Let us recall that we started our discussion from the fact that Gibbsian entropy, though 
intuitively acceptable, does not provide a measure of space-time order derived from a notion 
directly applicable to an individual configuration. At the same time, in the path-integral 
approach it is the order at the configuration level we are seeking when trying to identify the 
collective variables relevant in QCD. This led us to the information theory where the analog 
of Gibbsian entropy is represented by Shannon entropy, which explicitly relates to "ensem- 
ble" properties, and does not address the information content of individual objects from the 
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ensemble. 28 However, on the information theory side there exists a notion of Kolmogorov 
complexity which we used to define the "Kolmogorov entropy" of a configuration, and pro- 
posed that this is the appropriate measure of its space-time order. One can then define the 
corresponding measure for the Z 2 theory S via the ensemble average specified by Ps- We 
thus ended up with two distinct ensemble concepts for theory S, namely its Kolmogorov 
entropy and the Shannon entropy of the associated P$. What differences can one expect? 

It turns out that information theory provides a partial answer to this question for the- 
ories S that lead to computable probability distributions. The theorem bellow is a direct 
transcription of the relevant result for probability distributions on the set of binary strings 
(see e.g. Ref. [30]) to the framework developed here. 

Theorem 2 Let S G S GZ (for fixed size L of a symmetric torus) leads to a computable 
probability distribution P S (U), and E Sh [S] = -J2u p s(U) \nP s (U). Then 

< E K [S]-E Sh [S] < K(P S ) + 0(1) (32) 

where Ps(z(U)) = Ps(U) is the induced probability distribution on 23*. 

It needs to be emphasized here that the distribution P s (and its computability) has to be 
viewed as a function of the lattice size as well, i.e. Ps = Ps(%, L). We also recall that K(Ps) 
is the Kolmogorov complexity of P s (see discussion in section 3.3 and footnote 25). 

In view of the above theorem, we need to explicitly emphasize few points: 

(1) Theorem 2 implies that if Z 2 theories S\, S 2 have computable probability distributions, 
then their order relation induced by E K [S] is the same as the one induced by E Sh [S] m 
L — > 00 limit. Indeed, this follows from the equivalence relation (29) and the fact that 
E K [S] and E Sh [S] only differ by an L-independent constant for a theory with computable 
Ps- Thus, if it could be established that all theories in S GZ lead to computable distributions, 
one could also use Shannon entropy to distinguish different levels of space-time order in the 
infinite volume limit. In the specific case of Z 2 gauge theories this would justify the use of 
Gibbsian entropy (see footnote 28) as a measure of space-time order in the theory. However, 
the corresponding analysis of S GZ is not available yet, and one also needs to realize that 
the potential extension of such result would only be valid for theories with discrete local 
variables. 

(2) In connection with the above remark we wish to emphasize that we view E K [S] as the 
fundamental measure of space-time order in the theory, while E Sh \S] is an auxiliary concept 
that can be used, in some cases, to estimate E K [S\. For certain issues the use of E K [S] is 
simply insufficient. For example, consider the problem of minimizing E K [S] on some subset 
of S GZ with computable distributions. Unless K(P S ) is bounded on such subset (which is 
not necessary the case if the subset has infinite number of elements), then minimizing E Sh [S] 
will not help us at all even in the infinite volume limit. The recurring problem that appears 
in similar considerations is that Kolmogorov entropy and Shannon entropy are related by a 
"theory-dependent" constant (K(P S )). 

28 In fact, in case of Z2 gauge theories which we discuss, Shannon entropy is identical to Gibbs entropy. 
However, this is not true for probability distributions with continuous sample spaces (e.g. lattice QCD). 
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3.4 Kolmogorov Entropy: SU(N) Gauge Theories 

The extension of the above discussion to SU(N) gauge theories is straightforward in certain 
regards and problematic in others. From the point of view of numerical lattice QCD it is 
natural to view a configuration (to a given numerical accuracy) as a binary string since that's 
how it is in fact stored in the memory of a computer. While the basic idea remains to apply 
the Kolmogorov complexity measure to such strings and thus define the Kolmogorov entropy 
of a configuration/theory, the whole procedure needs to be formalized properly. Moreover, 
there are new conceptual problems arising due to the fact that local variables are continuous 
and thus involve infinitely many bits even on a finite lattice. The issues related to taking 
the appropriate limits are very relevant and need to be addressed. 

3.4.1 Canonical Coarse-Graining of the Gauge Group 

In this section we discuss a systematic way of discretizing the gauge group which is conve- 
nient for associating the SU(N) gauge configurations with binary strings. A particular way 
we propose is motivated by the issues of space-time order discussed here, but also by the 
anticipation of studying gauge theories from the point of view of information theory later. 
After discretizing the space-time coordinates via introduction of a lattice, it is quite natural 
to discretize all continuous aspects of the theory, including the continuous degrees of free- 
dom. However, this is usually not done since it is not necessary for proper regularization of a 
continuum system. At the same time, continuous gauge variables allow for incorporation of 
exact gauge invariance in essentially continuum-like manner. In reality though, exact gauge 
invariance carries mainly an aspect of "elegance" in this case since it is well under control 
even with a discretized gauge group, and its violations boil down to a numerical precision 
issue (See the discussion later in this section.). Moreover, for our purposes, the role of the 
gauge coarse-grained approximations is to define the measure of space-time order in the 
original lattice theory with continuous gauge variables, and not to supply a different lattice 
construction of continuum QCD. 

Discretization of gauge variables is made easier by the fact that SU(N) groups are compact 
and each element can be specified by N 2 -l real numbers chosen from bounded intervals. To 
provide a sufficient level of detail, let us focus on the group SU(3) which is most relevant for 
our purposes. It is convenient to start from a particular trigonometric parametrization for 
matrix elements of g e SU(3) in a fundamental representation, namely [32] 

gu = cos 6 1 cos 2 e i(f>1 g 13 = cos 0x sin 2 e^ 4 ^ 12 = sin^e^ 3 

g 22 = cos 0i cos # 3 e^ 2 g 32 = cos 6 1 sin 3 e 1 ^ 5 

021 = sin # 2 sin # 3 e^ 44 ^ - sin 0i cos 2 cos 3 e^ 1 ^ 2 "^ 

#33 = cos0 2 cos0 3 e- i( * 1+ * 2) - sin X sin 2 sin 3 e^ 3 ^ 4- ^ (33) 

g 23 = - cos 2 sin 3 e -^ 1+<fc) - sin 0i sin 2 cos 3 e^-^+W 

g 31 = -sin0 2 cos0 3 e- i(02+04) - sin0icos02sin0 3 e^ 1 ~ fe+ ^ ) 

< 0i, 2 , 3 < tt/2 < 0i,..., 05 < 2tt 
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with the invariant group measure given by 



dg = — - sin 9\ cos 3 O^dOi sin 2 cos 2 d0 2 sin 3 cos 3 d0 3 d<pi . . . d<p 5 (34) 

2"7T 5 

Using the above expressions we can now proceed to change the variables and define a canon- 
ical parametrization of the group where all eight parameters are chosen from a unit interval 
and the group measure becomes uniform (flat). In particular, we have 

Ui = 1 — cos 4 0\ u 2 = 1 — cos 2 2 u 3 = 1 — cos 2 3 u 3+a = — <p a (35) 

2ix 



where 



Since group elements are uniformly distributed in the space of parameters u a , we can 
proceed with the discretization in the usual manner. Instead of representing a value of u a 
exactly, we choose to keep only a "coarse-grained" information on it. In particular, to store 
such information in k bits, we can partition the associated unit interval into 2 k equal-sized 
subintervals and specify only which subinterval u a belongs to. Such arrangement obviously 
assumes the existence of a bijection between the set of subintervals and the set Q3£ of binary 
strings with length k. To set up such bijection let us first define a set V k of standard "lattice 
values" in the unit interval, and the "lattice spacing" e k for arbitrary k, namely 

V fe = {rf = ^, . = 0,1,2,. ..,2 fc -l} e k = 1 = rf> (37) 

Note that there is a natural bijection between *B* k and V k . Indeed, every i e N^ fe has a unique 
binary representation of the form % — b\ 2 k ~ 1 + 6 2 2 k ~ 2 + . . . + 2 + b k and then 

V k 3 rf } = 0.6x6a . . . b k <— b,b 2 ...b k = E *B* k (38) 

We will use the above map to define a bijection between binary strings and partitions of 
a unit interval by associating the values rf^ with individual subintervals. In doing so, 
we wish to treat periodic parameters («4, . . . ,u$) differently from non-periodic parameters 
(ui,U2,us) because of their different geometric meaning. In particular, we wish to end up 
with a discretization for which the immediate vicinity of an identity matrix corresponds to a 
single lattice cell rather than being split over several different ones. This requires the values 
of periodic parameters close to and 1 being associated with the same subinterval. We thus 
define 

23- = „«0 _ J I r i"'- r ?' + <*).* = 0, 1, 2* - 2 ,,„, 
" 3V - \ [rf'.rf + et ],i = 2--l {39 > 



for non-periodic parameters and 



« e k Jk) 



a, r ^ + f), * = 1,2,. ..,2* 



' [rf,f ]U[l-f,l), 2 = 
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for periodic ones. Let us denote the subintervals defined above as 3\ , with an implicit 
understanding of different definition for periodic and non-periodic cases. 

The bijections (38), (39,40) imply that we can interchangeably specify a value (r^), an 
interval (jf^) or a binary string {yf^) when speaking of "coarse-graining" a parameter u a 
into k bits. For future reference, let us define the operation of coarse-graining explicitly. In 
terms of a value of the parameter we have 

u a — ► u (k) = r (k) G V k such that u a G jf } (41) 

At the level of binary strings, we first point out that we can associate any possible value of 
u a with an infinite binary string z a = 616263 . . . formed by digits of its binary representation 
u a = O.616363 .... If u a has two binary representations (e.g. u a = 0.10000 . . . = 0.01111 . . .) 
we choose the option with infinitely many zeros. Note that we include u a — 1 via u a = 

0.11111 With these definitions the operation of coarse- graining into k bits translates at 

the level of binary strings into 

z a — 4*> ee yf > e n such that u a G t ] (42) 

Note that in the non-periodic case coarse-graining simply comes down to keeping first k 
bits of a binary representation for u a . In the periodic case it is the same after applying an 
appropriate periodic shift. 



3.4.2 Coarse— Grained Configurations and Binary Strings 

Consider an SU(3) configuration U = {U(n,fj,)} on a discretized hypercubic symmetric torus 
of size L. Emphasizing the canonical representation of the previous section, let us denote 
the same object as u = {u(n, [/,)}, where u(n, /i) ee (ui(n, //), U2(n, [/,), . . . , ug(n, //)). We 
associate with u (and hence with U) a sequence of fc-th level coarse-grained configurations 
u^ defined via 

u — u^ ee {«W(n,|i)} ee { (u[ k) (n, /x), «?> (n, /,),..., 4 fc) (n,/i))} (43) 

where the coarse-grained gauge parameters are defined in Eq. (41), and the u k-th level" refers 
to the fact that each parameter is represented by k bits. By construction lim^oo-u^ = u 
in the standard Euclidean norm. 

In an analogous manner we now proceed with defining a sequence of fc-th level binary 
strings corresponding to vS k \ To that effect we first define a binary string z^(n, fi) of length 
8k assigned to a given link (n, /i) via 

z^ (n, /i) ee z[ k) (n, zf (n, /*)... z^ (n, (44) 

where z a k \n, /i) are defined in (42) and the operation on the right-hand side means a simple 
concatenation. To form a string representing the whole configuration we have to fix the 
recursive invertible map F G T of Eqs. (23), relating a linear order on the string to underlying 
space-time geometry. For arbitrary such map we then have 

zf(XJ) ee z^\F-\\)) z {k \F'\2)) . . . 2 (fe) (F" 1 (4L 4 )) (45) 
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In other words, Zp is the concatenation of local binary strings z^ (n, /i) in the order provided 
by the map F. By construction, the length of this string is l(zp^) = 32L 4 k. 
Let us finally make few remarks pertaining to these issues. 

(1) Note that in addition to the choice of map F, there is additional freedom in assembling 
the binary string representing the coarse-grained configuration u^ k \ For example, instead 
of placing different bits of z^ (n, fx) into a localized part of the string, we could concatenate 
k blocks of size 32L 4 , where the j-ih block would be formed by the j-ih bit of z^ k \n,fx). 
However, these freedoms do not affect asymptotic considerations. 

(2) We emphasize that we will interchangeably use the notation and u^> to denote the 
same object. In the former case one should realize that U^(x,fi) is a matrix constructed 
from (x, n) via (33) rather than matrix of coarse-grained matrix elements. 

(3) For future reference, we will denote the set of all configurations U on the symmetric torus 
with L lattice sites on the side as U L . Similarly, the finite set of k-ih level coarse grained 
configurations will be denoted as U L,k C U L . 

3.4.3 Kolmogorov Entropy 

With the above preparations, we can now turn to the issue of quantifying the degree of 
space-time order for SU(3) (and analogously for SU(N)) configurations/theories. Following 
the same path as in the Z 2 case, the underlying idea is to apply the Kolmogorov complexity 
measure to binary strings describing the configurations. However, there are several problems 
with carrying out such a program. Among others, an immediate issue is that, even for finite 
L, the set X of objects we wish to describe (set of all SU(3) configurations) is not countable. 
This has several undesirable consequences, the most obvious one being that we cannot form 
a one-to-one correspondence between binary strings from *B* and configurations. 29 

In order to proceed, we will abandon the notion that the binary string associated with U 
has to encode this configuration completely. 30 Instead, we consider a sequence of coarse- 
grained configurations { U^ k \ k — 1, 2, . . .}, each of which approximates U in an increasingly 
fine detail. Clearly, the ability to analyze for arbitrary k is equivalent to analyzing U 
fully. At the same time, as we discussed in the previous section, each C/( fc ) can be completely 
described by a binary string Zp\lJ) of length 32L A k. For given fixed k, we can thus proceed 
and define the fc-th level Kolmogorov entropy of U corresponding to map F via 

E*(U,k) = K(zP(U)) (46) 

As discussed in Section 3.3, this measure is additively optimal with respect to set T with 
arbitrary F being an universal element. Following the same logic as in Z 2 case we fix F u of 
Eq. (26) and define 

E K {U,k) = K(z {k \U)) z {k \U) computed via F u (47) 

29 Note that it is not a problem to form a bijection between the sets of infinite strings and configurations 
at given fixed L. However, we need a bijection that includes all possible values of L. 

30 By "complete description" we mean such that any property of a configuration can be inferred from 
properties of the associated string and vice-versa. 
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If S = S(f3) is a theory from S G , we define its k-ih level Kolmogorov entropy as 

E K [S,k] = (E K (U,k)) s = J dUE K (U,k)P s (U) P s (U)=e' s ^/Z (48) 

i.e. as an ensemble average of Kolmogorov entropies assigned to individual configurations. 

In case of Z2 gauge theories, we could directly associate Kolmogorov entropy with the 
measure indicating the degree of space-time order dynamically generated by the theory in 
its typical configurations. For QCD the situation is slightly more involved. Indeed, it is 
important to realize that binary strings z^ k \U) can become asymptotically large by either 
increasing the lattice size L, or by increasing the level of coarse-graining k. Naively, one is 
tempted to take k —> 00 while keeping L fixed in order to completely characterize the theory 
in a finite lattice volume. While this can certainly be done, the Kolmogorov entropy in such 
limit will not be a measure of space-time order at all. Rather, it will reflect how ordered are 
the digits of a typical local gauge parameter u a (n,/i). Thus, it turns out that in order to 
obtain the construct we are interested in, we have to do just the opposite. In particular, we 
will say that E K (U, k, L) (E K [S, k, L\) is a measure of space-time order in the configuration 
(theory) if and only if 

k < L 4 (49) 

where we have made the L-dependence of E K explicit. 

The above considerations indicate that in dealing with issues of space-time order, we are 
simply required to approach any particular question at fixed coarse-graining level k. For 
example, the Kolmogorov entropy E K [S, k, L] induces a "Axth level" order relation on the 
set S G of theories defined on the torus with lattice size L, namely 

Si <k S 2 if E K [S 1 ,k,L] < E K [S 2 ,k,L] (50) 

The underlying interpretation is that Si possesses more space-time order than S 2 if Si <k S 2 
(for k <C L A ). Given that interactions are local, we expect that the above order relation 
becomes L-independent at sufficiently large L for any two theories, thus characterizing the 
interactions themselves. In the L — > 00 limit it is natural to consider an induced order 
relation on the set of equivalence classes, where theories whose Kolmogorov entropies differ 
at most by a constant belong to the same class (see Eq. (29)). 

Finally, let us point out that in complete analogy with the Z 2 case we can make the 
definition of k-th level Kolmogorov entropy gauge invariant by minimizing E K (U, k, L) on 
the gauge orbit of U (see Eq. (30)). We can also gauge-fix to the information gauge via an 
obvious modification of Eq. (31). 

3.4.4 The Conjecture of Hierarchic Structure 

One of the main reasons behind introducing the measure of space-time order based on 
Kolmogorov complexity is its universality in the sense that it can capture the ordered space- 
time behavior in any possible form. This is important because our geometric intuition for 
what the domain-correlated behavior of a function can be is very limited. Even taking just 
a scalar function of two space-time coordinates as an example that human brain can handle 
relatively comfortably visually, we have absolutely no intuition on possible degrees of space- 
time coherence beyond extreme cases. Indeed, we tend to deal with such functions via their 
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graphs which we can picture mentally if they are nice smooth surfaces in 3-d space. The 
common intuition then is that smoother graphs (functions) are "more ordered" than rougher 
graphs (functions). However, even in the extreme case of very high order our intuition is 
incomplete since there exist very singular functions whose behavior is maximally ordered in 
the sense that local patterns determine their behavior everywhere. 

Returning back to lattice configurations, the above-mentioned generality of the approach 
based on Kolmogorov complexity, while highly desirable, has certain implications that have 
to be dealt with. In particular, the point that needs to be emphasized at this stage is 
that the linear order on the set of configurations induced by E K (U, k, L) is not necessarily 
preserved if the "resolution" k changes. More precisely, if U\, U 2 G U L , and k\ < k 2 <C L A 
with U\ <fci U 2 , then it does not necessarily follow that U\ <k 2 U 2 . Indeed, we can easily 
imagine a situation where first k\ bits of gauge parameters in U\ exhibit high space-time 
coherence, while being disordered in U 2 . However, in the trailing k 2 — k\ bits the situation 
can reverse in favor of U 2 if supplied with highly ordered digits. The point is that while 
we intuitively tend to think that leading digits are "more important" for the structure, the 
Kolmogorov complexity measure does not distinguish such hierarchy and the order around 
1-st digit counts in principle as much as the order around the 100-th digit. We emphasize 
again that this is a desirable feature of Kolmogorov entropy measure since our goal is to 
be completely general. Thus, if there is a secret message encoded into configuration U via 
inserting a highly space-time ordered binary string into various digits of u, it can still be 
detected by analyzing E K (U, k, L) for sufficiently large k. 

The point we wish to convey in the above reasoning is that at the configuration level, 
the ranking introduced by Kolmogorov complexity has to be viewed strictly as a resolution- 
dependent concept in this case. In other words, it is meaningful to say that configuration U\ 
is more ordered than configuration U 2 only when they are both analyzed with fixed resolution 
k. However, when we characterize the theory S(g) G S G , via E K [S, k] things are expected 
to be much more universal in this regard. Indeed, since our theories are analytic (at least 
in certain corner of the configuration space), we expect that with other things being equal 
(such as multiplicities), the space-time order in leading digits has a more decisive influence 
on probability of occurrence than possibly equivalent order in higher digits. This is simply 
because leading digits influence the value of the action (and thus a probability of occurrence) 
more significantly than higher digits. In other words, it is natural to assume that a theory 
in S G will feature typical configurations with hierarchical order. By this we mean that 
the relative degree of order will not increase when higher digits are included, and will in 
fact typically decrease. Put still differently, the higher digits in typical configurations are 
expected to exhibit increasingly more space-time randomness. 

We wish to formalize the above observations into a definite statement concerning the 
measure E K [S,k] of space-time order on the set S G . While at present we do not have 
an apparatus to attempt the proof of the following statement, we propose it for future 
considerations because of its relevance. 

Conjecture CI1. Let k be a fixed positive integer and let "<k" be the linear order relation 
on S G induced by E K [S, k, L] in the L — > oo limit. Then "<k" is universal with respect to 
k, i.e. if Si <ki S 2 , then S± <k 2 S 2 , for arbitrary ki, k 2 and arbitrary Si, S 2 G S G . 
Few points need to be emphasized in connection with the above conjecture. 
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(1) The utility of Kolmogorov entropy as a measure quantifying the degree of space-time 
order in different theories does not rely on the validity of CI1. Indeed, in our quest for 
uncovering the relevant collective degrees of freedom we can use E K [S, k] at fixed resolution 
k as a general measure distinguishing different theories. Obtaining a complete picture will 
require changing the resolution k in either case. Nevertheless, from a theoretical standpoint it 
would be very appealing if CI1 could be demonstrated at least in some simplified framework. 
In what follows we will implicitly assume its validity. 

(2) A weaker statement with similar consequences (namely existence of an absolute ranking 
of theories by the degree of space-time order they generate) can be obtained if one postulates 
the existence of k Q [Si,S 2 } such that the relative ordering of Si and S 2 is preserved for all 
k > k . 

(3) We remind the reader again that it is crucial in the formulation of CI1 that the infinite- 
volume limit is taken at fixed k to define the ordering. This is to nullify the influence of 
various non- universal constants implicitly present in the definition of E K [S, k, L\. Also, it 
is implicitly understood in CI1 that in L — > oo limit the linear order is defined on the 
equivalence classes of theories whose Kolmogorov entropies are bounded by a constant (see 
Eq. (29)). 

3.4.5 The Associated Coarse-Grained Theory and Shannon Entropy 

In section 3.3.2 we have argued that in the case of Z 2 gauge theories it is useful to consider 
the Shannon entropy of the corresponding ensembles. Even though this is an auxiliary con- 
cept from our point of view, its usefulness stems from the fact that Kolmogorov entropy is a 
non-computable notion. Using Theorem 2, Shannon entropy gives us at least a possibility to 
calculate which of the arbitrary two Z 2 gauge theories with computable probability distri- 
butions generates more space-time order in the infinite-volume limit. This raises a natural 
question at this point, namely how to define this auxiliary concept in case of SU(N) lattice 
gauge theories. 

As we argued extensively above, we carry out the analysis of space-time structure at 
a given coarse-grained level k, wherein the configuration U is associated with its coarse- 
grained counterpart (or u'*'), completely represented by a binary string z^ k \ Given 
S E S G and fixed L, our goal is to define a probability distribution on a finite set U L ' k of 
coarse-grained configurations U^ k \ such that its Shannon entropy is related to Kolmogorov 
entropy E K [S,k,L] in the same manner as in Theorem 2 (see Eq. (32)). To do that, it is 
useful to think about the theory S in terms of the corresponding ensemble S[S], which is 
fully represented by an appropriately chosen infinite probabilistic chain. The coarse-graining 
procedure corresponds to the replacement 

£[S] = {...,[/,_!, £/;,£W..} — £^[S] = {...M k lM k) M%...} (51) 

The ensemble defines a lattice theory of discrete variables which we denote 

as S( k \ so that we can write formally £^[S] = £[S^]. 31 The probability density of 

31 Dcfinition of action is more than formal though in the sense that its value for arbitrary £ U L,k 
can be explicitly written down (see bellow). 
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encountering a configuration U G IA L in £[S\ is given by Ps{U) = exp(—S(U))/Z , while the 
probability of encountering [/^G W L ' fc in £^[5'] is by construction 

P 5(fc) ([/W) = J dV5(u {k) -V {k) )P s (v) (52) 

where is a fc-th level coarse-grained configuration associated with V. Note that using 
our canonical representation of gauge variables, we can write the above equation in a precise 
(rather than formal) manner as 

P SW (U^) = [ n dv a (x,n) J] 5(ui k \x^)-v^(x^))P s (v) (53) 

Moreover, for actual evaluation of this probability one could use an explicit form without the 
5-functions. To do that, we define the k-th level displacement variables = {t^(x, pi) } 
with values in the intervals 

r {k) (x u)el '°' e ^' a=l,2,3 , 

r a [_ ejfc/2> Cfc / 2 ] , a =4,...,8 

in terms of which we can write 

where the periodicity is implicitly assumed to be taken into account for arguments of Pg 
with a = 4, . . . , 8. Equation (55) defines a theory of discrete variables such that its 
Kolmogorov entropy is equal to k-th level Kolmogorov entropy of S", namely 

E K [SW\ = Y, EK ( u{k) ) p sw(U {k) ) = E K [S,k] = I ' dUE K (U,k)P s (U) (56) 

where in the context of (left hand side; discrete variables) we have E K (U^) = K(z^), 
which is equal to E K (U, k) in the context of S (right hand side; continuous variables). 

For theory with computable probability distribution, there is a relation between its 
Kolmogorov and Shannon entropies analogous to Theorem 2 in Z 2 case. This then induces 
the same relation between Shannon entropy of and the fc-th level Kolmogorov entropy 
of S which is summarized in the following theorem. 

Theorem 3 For L, k fixed, let S E S G leads to a discrete theory with computable 
probability distribution P s(k) (U^), and E Sh [S^] = - £i/W P s w(U^) lnP 5(fc) ([/«). Then 

< E K [S,k] -E Sh [S (k) ] < K(P SW ) + 0(1) (57) 

where P S (k){z = z^) = P S (k)(U^) is the induced probability distribution on *B*. 
To close this section, we wish to emphasize two points. 

(1) The remarks given at the end of section 3.3.2 apply in the case discussed here as well. 

(2) The associated coarse-grained theory is not exactly gauge-invariant with violations 
at the k-th bit level. This does not have any undesirable implications in this context since 
the role of is to assist us with quantifying the degree of space-time order in S. 
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4 The Principle of Chiral Ordering 



In the previous section we have introduced the notion of Kolmogorov entropy for SU(N) gauge 
theories and proposed it as a measure quantifying the degree of space-time order generated by 
the theory at the configuration level. Since the concept is based on Kolmogorov complexity 
of binary strings describing coarse-grained configurations, it is completely general in the 
sense that it reflects space-time order of coarse-grained configurations/theories in any form. 
Moreover, if we accept the hypothesis of hierarchic order (see conjecture CI1), then we have 
a measure providing universal ranking of theories in S G by the degree of space-time order 
generated in their typical configurations. 

Before proceeding further, let us note that we included a detailed elaboration on Kol- 
mogorov entropy in the previous section simply because we view the recent findings expressed 
in II (see section 1.3) as an important new aspect that needs to be properly incorporated into 
the framework of quantum field theory. Equipped with the notion of Kolmogorov entropy 
we can, at least at the conceptual level, attempt to formulate problems that are relevant for 
the program of uncovering the fundamental structure of QCD vacuum in the path integral 
formalism. In that regard, it we would certainly be of crucial importance if we had means 
to achieve the following major goal. 

Goal 1 : Map out the landscape of Kolmogorov entropy in the set of actions S G and under- 
stand its behavior. 

Indeed, in order to uncover the nature of relevant collective variables C, we would eventually 
prefer to deal with theories possessing low Kolmogorov entropy at finite lattice spacing. In 
such theories the distorting influence of randomness on C is expected to be low thus giving 
us best chance of correctly identifying it. 

To start speaking in these terms, let us clear up some terminology first. While S G is 
the set of individual lattice actions, we tend to speak of lattice "theories" belonging there 
in two different ways. Firstly, we view the quantum dynamics of finite (countable) degrees 
of freedom defined by any individual element S(g) (g and other parameters fixed) of S G 
as a "theory". Note that we talked about it in this sense when discussing the Kolmogorov 
entropy. For example, we can compare the space-time order of two theories differing only 
by values of g, e.g. S(gi) <k S(g 2 ). 32 On the other hand, it is sometimes useful to think of 
lattice theory as a one-parameter family of interactions S(g) parametrized by bare coupling. 
In this case the "theory" prescribes a lattice interaction at any value of g thus parametrizing 
the approach to the continuum limit. 

Thinking now in terms of this second interpretation, consider the set of pure-glue lattice 
gauge theories contained in S G . Upon quantization, each theory S(g) G S G is associated 
with ensembles £(g) which determine the lattice quantum averages for various observables. 
To obtain values in physical units, one usually fixes the appropriate quantity to its physical 
value, independently of g, which introduces the lattice spacing a. From the physical point 
of view it makes clearly more sense to compare space-time order of two theories at the 
same lattice spacing rather than the same bare coupling, i.e. to view the theory S as being 
parametrized by lattice spacing rather than bare coupling. Since this is what we will do, let 

32 In fact, we naturally expect that typical configurations will become increasingly spatially ordered as g 
decreases. In other words, we expect that S(gi) <fe S(g2) for any gi < g<2. 
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us describe the procedure of introducing the lattice spacing in more detail. Assume, for the 
sake of argument, that the mass of the lightest scalar glueball has been fixed to its value rh 
for all theories S(g), in each of which its dimensionless value is m s (g). For given S(g) this 
defines the lattice spacing a s (g) via 

a s {g)m = m s {g) (58) 

For simplicity, assume that the set of other interesting physical observables consists of masses 
rhi, of all the other stable particles in the theory, which we want to predict. In this restricted 
setting the statement that all actions S(g) G S G define the same continuum theory (univer- 
sality) means that 

lim — I, . = m lim — l -^-rr = Vi, independent of S (59) 

9^0 a s (g) Hm s (j) 

In other words, all lattice theories S G S G will predict the same ratios of masses of stable 
particles in the continuum (g — > 0) limit. Note that in the above considerations we have 
implicitly assumed that all theories in S G have a second order phase transition at g — 0. If 
we further assume that the lattice correlation lengths £f(<?) = l/mf(g) are monotonically 
increasing as g — > + (at least for g < g, where g is independent of S), then we can invert the 
functions a (g) and talk about lattice actions parametrized by lattice spacing a. In other 
words, we set S'(a) = S(g = g s (a)), but in what follows we will keep the same symbol (S) 
for both dependences to simplify the notation. 

With all the above pieces in place, it is tempting to conclude that, at least at the concep- 
tual level, it is now clear how to proceed with the program of identifying the fundamental 
structure of QCD vacuum (collective variables C). Being restricted to some finite range of 
non-zero lattice spacings a by practical considerations, we could minimize the Kolmogorov 
entropy E K [S(a), k] in the set of actions at fixed a to identify the action S°(a) which max- 
imizes the space-time order. We would then study the corresponding ensemble S°(a) in a 
Bottom-Up manner to extract as much information as possible about the nature of space- 
time order present. The lattice spacing could then be decreased and the procedure repeated 
in an attempt to gain the information about the sought-for collective variables C in the 
continuum limit. 

As simple and appealing this might sound, it is most likely a wrong way to pose the 
relevant task (even at the level of principle, leaving out the practical considerations). Indeed, 
apart from the issues of existence and uniqueness of S°(a) (which might not be serious in 
this case 33 ), there is a fatal flaw in the program as formulated above. Indeed, assume that 
we start the minimization to determine S°(a) from the action S im (a) such that theory S mi 
is well in the scaling region for lattice spacing a considered. Thus, for S mi we have that 
mf m \a)/a fa rhi for physical masses m t . Since the minimization of E K makes no reference 
to correlation lengths other than £ 5 (a) = l/m s (a) = 1/arn (recall that lattice spacing is 

33 We note that such considerations were actually part of the reason why we chose to parametrize theories 
in terms of a rather than g. Indeed, in the latter case, the set S G restricted to actions at fixed g most likely 
does not have a bottom element with respect to E K . The other reason obviously is that when we compare 
Kolmogorov entropies for two actions with the same number of lattice degrees of freedom and fixed a, we 
are comparing the degrees of space-time order in the same physical volume. 
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kept fixed), there is no reason to expect that the scaling behavior will be preserved during 
the minimization. In other words, the masses of stable particles (other than the one fixing 
the lattice spacing) would most likely drift away from their physical values. Thus, even 
though the Kolmogorov entropy might be low in S°(a), we would certainly not want to use 
it to make conclusions about the nature of space-time order in the continuum limit, no 
matter how small a might be. These considerations suggest that in an attempt to develop 
a systematic program for studying space-time structure in QCD vacuum, we would ideally 
(in principle) like to be able to do the following. 

Goal 2 : Starting from some action S mi (a), we would like to move in the set S G toward 
actions with lower Kolmogorov entropy such that if S im {a) is in the scaling region, then the 
theories reached will remain in the scaling region, and if S mt (a) is not in the scaling region, 
then the theories reached will not diverge further away from it. 

As we emphasized along the way, the above goal specifies the "in principle" strategy for 
identifying a fundamental structure in QCD vacuum. Indeed, it appears unlikely that at the 
practical level we could proceed via direct constrained minimization of Kolmogorov entropy 
in the set of actions. Instead, we will view Goal 2 as a conceptual guide against which our 
attempts should measure their potential merit. The second part of this article deals with a 
proposal in that direction. 

To start describing our approach, it is instructive to first highlight the crucial obstacles 
to practical realization of Goal 2 which need to be addressed, (a) There is no standard 
explicit way of parametrizing the elements of S G and thus no standard practical way of 
"moving in the set of actions" . There is no clear underlying principle that would guide 
us to move in S G while "preserving physics" as demanded by Goal 2. One could view the 
renormalization group (RG) blocking as a procedure of such kind, except that the decimation 
of the lattice poses a problem. At the heart of our method is the proposal for such principle. 
(7) Since Kolmogorov entropy is not computable, we cannot directly minimize it or at least 
check aposteriori by explicit calculation that it became lower during the procedure. One 
thus needs to discuss why is the chosen evolution in the set of actions believed to lower the 
Kolmogorov entropy. In the sections that follow we will focus on the above issues in turn. 

4.1 Configuration-Based Deformation of the Action 

In this section we describe a generic way of "moving" in the set of actions based on the local 
deformation of individual configurations. We will use the fact that it is not important for us 
to know analytically (or otherwise) the exact form of the action we work with. Indeed, we 
only need to be able to inspect the ensemble associated with the action and to make sure 
that it represents a valid element of S G . Thus, consider a map 

•M. : U L — ► U L U r ^ — > M n>li (U) (60) 

Let S G S G and Ss be the associated ensemble. The map Ai then induces the new ensemble 

S s = {...U l '\U\U %+ \...} — > {...7W(^- 1 ),7W(f/ i ),7W(f/ m ),...} = £ s m (61) 

where the old probability density Ps(U) = e~ s ^ jZ$ is modified to 

P SM (U) = j dV P S {V) 8(U - M(V)) (62) 
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and is associated with the "action" S M given by 

e -sM(u) = p sM (jj)Z S M Z s m = Z s (63) 

We now impose the following requirements on the map M.: 

(i) Symmetries. We require that the symmetries of the original action S(U) G S G are 
also present in new action S M (U). One can easily inspect that this is guaranteed if M. n ^ 
transforms in the same way as U n ^ under respective symmetry transformations. Thus, for 
example 

M n ,,(U) -^G n M n ,^U)G- 1 +[i under U n>li — G n U n ,,G~^ (64) 
and similarly for the symmetries of the hypercubic lattice structure. 

(ii) Locality. Considering the configuration in the canonical representation (see Sec. 3.4.1), 
i.e. U = {Un tfl } = {un tfl } = u we require the existence of constants C, 7 such that 

dM Y {u) < Ce-^ n - m \ (65) 

independently of u. 

(Hi) Classical Limit. If A^x) is a classical (smooth) field (see footnote 9), and U n>fi (a) is its 
discretization (9) with classical lattice spacing a, then we require that M. n>fi (U, a) is also its 
valid discretization. In other words, that 

U n M-M n ^(U,a) = 0(a 2 ) (66) 

Note that the nature of the above conditions (especially (ii) and (Hi)) is such that it is 
implicitly assumed that M. is defined for arbitrary L. In the following, we will denote the 
set of maps satisfying the conditions (i)-(iii) as 3". 

We propose that the action S defined via Ai G 3" satisfying (i)-(iii) will generically be 
a valid element of S G . It needs to be emphasized that this conclusion is not the result of 
a rigorous proof at this stage. Indeed, while the symmetries are directly guaranteed by (i), 
it is not straightforward to precisely formulate and prove the relevant statement regarding 
locality of S based on (65), or the statement on the classical limit based on (66). 34 On 
the other hand, from the physical point of view it is quite clear that the transition from S to 
S M will not cause any major pathologies if (i)-(iii) are satisfied. Indeed, consider arbitrary 
set of local operators {O a } and the related set {O^ 1 }, where 

0»(U) = O a (M(U)) (67) 

The conditions (i)-(iii) guarantee that will inherit the transformation properties, locality 
and classical limit of O a , while at the same time 

(O a (U)) 8M = {0?(U)) S (68) 

Thus the expectation values of arbitrary local operators (and their correlation functions) 
in S correspond to values in original theory S (element of S G ) for manifestly "valid" 
discretizations of the same continuum operators. Roughly speaking, our overall level of rigor 
here is analogous to that of Wilson's renormalization group (RG) approach [33] , wherein one 
moves in the set of actions via RG transformation. 



34 This is not to say that it is impossible to formulate and prove such statements (with appropriate restric- 
tive conditions stipulated). However, we will not deal with these rather involved issues rigorously here. 



36 



4.2 The Principle of Chiral Ordering 

Consider a configuration-based deformation S M of action S, given by an ultralocal map 
M. G £F. By ultralocality we mean that M n ^(U) does not depend on XJ m ,v for |n — m\ larger 
than some configuration-independent positive constant k. Using the same arguments as 
those invoked in the case of RG transformations, we can conclude that the behavior of S 
at large lattice distances is the same as that of S. In other words, the bound-state masses 
rhf in lattice units are expected to be the same (on an infinite lattice) as fhf M . However, 
there are several setbacks associated with such generic approach to moving in the set of 
actions, (a) Transition from S to S M can strongly affect observables for which the gauge 
field fluctuations at the scale of k lattice spacings are crucial (i.e. observables that are not 
determined entirely by long-distance behavior of correlators), even on an infinite lattice. 
(b) On a finite lattice, even the estimates of masses in S M can change from the estimates 
in S. Moreover, after applying the deformation many times, all observables in S M " on a 
finite lattice can change uncontrollably (both in lattice and physical units) relative to S. 
There is nothing to protect them from changing, (c) For the generic choice of Ai, (which 
typically can be interpreted as a "smoothing" operation), the repeated application of the 
map on given configuration (A4 k (U)) is not expected to converge to a non-trivial limit. 
Indeed, one rather expects the convergence to a trivial configuration with gauge potentials 
vanishing (U n ^ = I (identity matrix)), or a diverging behavior. Thus, even if Kolmogorov 
entropy of M.(U) decreases relative to U for all relevant configurations, (which also means 
that E K [S M ,L] < E K [S, L]), we cannot expect a unique/interesting result of minimizing 
procedure based on the repeated application of M.. Taken together, the above points serve 
to illustrate that to approach Goal 2 via configuration-based deformation of the action with 
generic M. is not expected to succeed. We also note that these issues become slightly more 
complicated in the non-ultralocal case since lattice correlation lengths in S can in principle 
change (by at most the range of non-ultralocality) relative to S, and thus the masses (and 
the lattice spacing) can change accordingly even on an infinite lattice. 

In this section we wish to propose a class of maps M. which, we believe, can potentially 
alleviate most of the problems encountered in generic case. It is quite obvious from the 
above discussion that such maps must be a very special subset of all possibilities and, in an 
attempt to identify them, we will proceed quite heuristically. Indeed, our only initial guide 
will be the striking success of linking the topological density operator to chiral symmetry, 
which resulted in II of Sec. 1.3, i.e. the observation of space-time order at the fundamental 
level [5]. It will become clear soon in the process however, that the approach potentially 
carries much deeper connotations. Indeed, we will use chiral symmetry here in different and 
more general manner. Consider a fermionic action S F specified by lattice Dirac kernel D. 
Apart from all conditions imposed on D(U) by the fact that S F is an element of S F (which 
includes requirement of chiral symmetry), we impose an additional and highly non-trivial 
requirement. In particular, 

D n , n+ ,(U) = D[ x Vn, n , VC/ U n „ G SU(N) (69) 

where is a free fermion matrix (no color indices) representing hopping by vector m, 
namely 

{D f m ) a ^ = D ntn+m (l) Q)/3 color index a fixed. (70) 

a, a 
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In other words we require that the color-spin matrix D ntn+tl (U) is a direct product of a 
free fermion matrix D* and an SU(N) "matrix phase" U ntfi . Consequently, the fermionic 
action S F specified by D satisfying (69) defines a map M. D which associates with every 
configuration U a new configuration U via 

U — > M D {U) M^(U) ee U n ,» (71) 

We will call U = M. D {U) a chirally ordered version of U according to D. One can easily 
check that the map M D satisfies the conditions (i)-(iii) of Sec. 4.1. 35 Consequently, M D 
is expected to induce a valid deformation of the gauge action S — ► S M G S G . Action 
S M will be referred to as a chirally ordered version of S according to D. 

Few crucial points pertaining to the above definition of chiral ordering need to be made clear 
at this point. 

(1) While it is not straightforward to prove that the elements of S F satisfying (69) exist, 
their existence does not appear to contradict any requirements imposed on valid fermionic 
actions. 36 Moreover, later in this paper, we will claim that there exist converging sequences 
of elements in S F (with the limit that could be outside the set) that will satisfy (69) to 
arbitrary precision. Accepting conditionally the possibility, let us denote the subset of such 
special fermionic actions as S F ' C . Note also that from now on we will frequently identify the 
action S F with its defining kernel D and so we will interchangeably write e.g. S F G S F and 
D G S F . 

(2) According to Eq. (69) the operation of chiral ordering replaces the link U n ^ with the 
interaction matrix phase U n ,n involved in "hopping" of chiral fermion (quark) from site 
n + \i to site n. In other words, it is a measure of how much is the elementary hopping of 
chiral fermion affected by the presence of non-zero gauge field relative to the dynamics of free 
fermion. Thus, as far as chiral fermion is concerned, the physics content of configuration U is 
very much embodied in configuration U = M. D {U). This is our rationale for the expectation 
that the relevant physics encoded in gauge action S is very close to physics in S M . We 
take this as a basis for the proposition that moving in the set of actions via chiral ordering 
can potentially solve the problem {(3) stipulated in the introduction to Sec. 4. 37 

(3) We emphasize that the condition (69) can be easily satisfied by lattice fermionic actions 
that are not chirally symmetric. For example, the Wilson-Dirac operator Dw (which is not 
an element of S F ) satisfies it with U n ^ = U n ^, i.e. with M Dw being an identity map. It 
is thus the chiral symmetry of the probing fermion that turns the chiral ordering condition 
into a highly non-trivial requirement. This is natural since we view the chirality of lattice 
fermion as a necessary requirement for properly regularizing the continuum physics. At the 
technical level, this comes as a result of non-ultralocality [34] of D G S F in gauge variables. 
Indeed, as a consequence of Ginsparg-Wilson relation, the matrix D ntn+fl (U) is expected 
to receive contributions from all possible gauge paths starting at n and ending at n + /i. 

35 Indeed, the locality and the symmetries follow from the locality and the symmetries of S F . The classical 
limit is also preserved as one can see from the expansion of D n ^ n+fl (U) in terms of gauge-covariant paths. 

36 The non-trivial part in constructing such elements obviously comes from reconciling the requirement of 
chiral symmetry (Ginsparg-Wilson relation) with equation (69). See also point (3). 

37 See also points (a)-(c) in the beginning of this section. 
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Condition (69) demands that these contributions add up to an effective SU(N) phase. This 
is clearly a highly complex requirement, and the associated M. D is a highly complex map. 

(4) Viewed as a requirement for the fermionic action, Eq. (69) represents a condition that 
is independent from the usual ones imposed on elements of S F . In particular, it is not 
implied by standard symmetries. 38 This is why we sometimes use the term "principle" 
when referring to it. However, our present aim in imposing Eq. (69) is not to define a better 
lattice fermionic theory, but to provide a guide for evolution in the set of gauge actions that 
suits our goal of studying the QCD vacuum structure. In that regard, we view the principle 
of chiral ordering as the following generic statement. 

The Principle of Chiral Ordering: The configuration-based deformation of the gauge 
action defined by map M. representing the effective SU(N) interaction phase associated with 
elementary hopping of lattice chiral fermion has a very mild effect on lattice physics. This 
applies to both long and short-distance behavior (in lattice units). 

Thus, what we mean by the principle of chiral ordering is the broad statement according 
to which using the maps Ai n ,fi with tendency to capture the effective interaction phase of 
chiral fermion are particularly suitable for exploring the nature of space-time order in QCD. 
We emphasize that the potential usefulness of this observation does not crucially depend 
on the fact that the set S F ' C C S F is not empty. Indeed, we view the condition (69) as a 
clearest expression for the underlying idea, which we can use as a guide for extracting the 
approximate matrix phase corresponding to arbitrary chiral lattice fermion. This will be 
discussed more fully in Sec. 5. 

4.3 The "Ordering" Conjecture 

In the previous two subsections we have addressed the points (a), {(3) put forward in the 
opening part of Sec. 4. In particular, we have suggested how to move in the set of actions (via 
configuration-based deformation) while minimally changing the lattice physics (via principle 
of chiral ordering). In these considerations we have already implicitly assumed, as the name 
"chiral ordering" suggests, that this will also lead to lowering of the associated Kolmogorov 
entropy, i.e. that the point (7) is automatically satisfied. To prove the relevant rigorous 
statement within the framework of Kolmogorov complexity is a difficult task that we will 
not attempt here. However, we need to make it clear that our considerations assume the 
validity of the following conjecture. 

Conjecture CI2. Let k be an arbitrary but fixed positive integer. Then E K [S M ° , k, L] < 
E K [S, k, L] for sufficiently large L. 

In other words, we assert that chiral ordering leads to a gauge theory with lower Kolmogorov 
entropy. Note that the only restriction on lattice size L in CI2 is related to applicability 
of Kolmogorov complexity concepts. Also, we do not restrict the validity of CI2 to chiral 
orderings defined by D e S F,C , i.e. we propose that it is valid also for the approximate 
orderings performed with arbitrary D G S F , which we will discuss later in this paper. It 
should also be emphasized that we do not claim that the conjecture is strictly valid at the 

38 This does not necessarily mean that the condition cannot be phrased in the language of symmetry. 
However, if it does, then the associated symmetry is different from the standard ones. 
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configuration level, i.e. it is not necessary true that E K (M D (U), k) < E K (U, k) for arbitrary 
U G U L , even though it is reasonable to expect that it is true for most configurations. 

The heuristic argument for validity of the "ordering" conjecture is tied to the principle 
of chiral ordering and is as follows. The underlying idea behind the principle of chiral 
ordering is that the physical content of a configuration mostly resides in the SU(N) interaction 
phases associated with elementary hopping of a lattice chiral fermion. At the same time, 
the path integral approach to QCD vacuum implicitly relies on the assumption that the 
fluctuations associated with "physics" exhibit high degree of space-time order (and will lead 
to identification of collective variables C), while the physically irrelevant fluctuations are 
mostly disordered. Thus, by identifying and keeping the physically relevant part of the 
fluctuations, we should increase the order on average, and thus to lower the Kolmogorov 
entropy. Needless to say, it would be very illuminating if the statement analogous to CI2 
could be proved rigorously at least in some model situation. 

4.4 Repeated Chiral Ordering and Perfectly Chirally Ordered 
Configurations 

It is interesting (and in the case of approximate chiral ordering very important) to think 
about repeating the chiral ordering procedure iteratively. To start speaking in these terms, 
let us denote the configuration obtained from U in j-th step as More precisely 

jjio) _ u _ ( M D ) j (U) , j = 1,2,... (72) 

At the same time, the original gauge action S = S^ D evolves into after j steps. 

The configuration U^(— U) represents the interaction-phase content of U^\— U). Is the 
interaction-phase content of equal to itself, i.e. is equal to U^? Even though this 
would seem natural, it is not necessarily so. On the other hand, from the point of view of the 
principle of chiral ordering, the configurations stable under map Ai D are clearly special. In 
fact, we expect them to be the configurations that reflect the space-time order present in the 
vacuum in the clearest possible manner. Let us thus define the corresponding set precisely. 

Definition 1. Let D{U) G S F be an arbitrary chiral lattice Dirac operator defined on a 
latticized torus of size L (i.e. U G U L ). The configuration U satisfying the condition 

D n , n+ ^(U) =Dfx Un^ Vn, » (73) 

is said to be perfectly chirally ordered. The subset of all such configurations will be denoted 
as U L ' D . 

We emphasize that in the above definition we do not restrict ourselves to Dirac operators that 
are elements of S F ' . Indeed, the set of perfectly chirally ordered configurations is defined 
for arbitrary chirally symmetric action, and will play an important role in the procedure of 
approximate chiral ordering discussed later in this paper. 

As is obvious from its definition, the set U L,D is not empty for arbitrary D, since it 
contains a trivial configuration with U n ^ = I c ( "free-field" configuration; I c is an identity in 
color space). However, the notion of perfectly chirally ordered configuration (and the whole 
program of chiral ordering) would not be of much interest if this was the only possibility. 
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Thus, the following conjecture is to be viewed as one of the cornerstones for the approach 
to vacuum structure proposed in this paper. 

Conjecture CI. The set U L,D contains non-trivial configurations for arbitrary D G S F 
for sufficiently large L. More precisely, it contains configurations with non-vanishing gauge 
invariant composite fields. 

An important feature of the Conjecture CI is that one can attempt its numerical verification 
for overlap fermions [12]. 

Returning now to the repeated exact chiral ordering, let us consider the sequences 
{U®,j = 0,1,2,...} and {S$ D ,j = 0,1,2,...} defined above for D e S F ' C . Accord- 
ing to the principle of chiral ordering and Conjecture CI2 we expect that, with increasing j, 
we increasingly reduce the randomness present in the configuration, while changing lattice 
physics very mildly. Thus, in some sense, our approach expresses the view that random 
fluctuations are mostly responsible for violating the condition (69), while the relevant fluctu- 
ations ( "physics" ) tend to satisfy it. For the elements of U L ' D this process ends immediately - 
they are invariant configurations under M. D . For U ^ U L,D it is expected that the amount of 
removed noise decreases fast with j and, asymptotically, it is not possible to remove any more 
randomness. In other words, the elements of { U^\j = 0, 1, 2, ... } behave asymptotically 
as elements of U L,D . This expectation is expressed in the following conjecture. 

Conjecture C2. Let U ElA L and D e S F,C . (a) The limit lim^oo exists generically 
and is an element ofU L,D . (b) There exist configurations U for which this limit is non-trivial 
at sufficiently large L. 

Note that by "generically" we mean up to some possible small set of measure zero which 
we postulate not to be relevant in the continuum limit. We wish to emphasize the following 
points. 

(1) The non-trivial part of conjecture C2 is the generic existence of the limit and the fact 
that there exist configurations for which the limit is non-trivial. Indeed, if lim^oo = U 
exists, then it is easy to see that it must be an element of U L ' D since 

lim || U® - MZjpV) || = lim || U® - \\ = 

and so Ai^^U) = U n ,n- In the above equations ||.|| denotes an arbitrary matrix norm. The 
existence of a non-trivial limit (at the configuration level) for the procedure of chiral ordering 
is one of the main features that, we believe, distinguishes this approach from a deformation 
based on a generic map M. (See point (c) in Sec. 4.2.). In a way we are suggesting that the 
deformation that preserves the interaction phases (and hence the physics) avoids converging 
to trivial configurations where physics is clearly changed. Instead, the tendency to preserve 
physics is strong enough so that a convergence to non-trivial invariant configurations arises. 

(2) As & result of conjecture C2 the sequence { S^ D ,j = 0, 1, 2, ... } converges to a non- 
trivial action S(S, M D ). By ' 'non-trivial" we mean that the corresponding action assigns 
a non-zero probability density to at least some non-trivial configurations. Note that while 
all individual actions are expected to be valid elements of S G , we do not necessarily 

claim that the limiting action S is also such an element. Indeed, while the classical limit and 
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the symmetries are surely preserved, the action could in principle become asymptotically 
non-local. However, we emphasize that even if this happened in some cases, the relevant 
message is that S can be approximated to arbitrary precision by an element of S . 

(3) At the operational level, the special significance of the elements of U L,D (invariant con- 
figurations) comes from the fact that their probability of occurrence for actions reached via 
repeated chiral ordering cannot decrease. In fact, the total probability of these configura- 
tions and their small deformations will grow and, according to the conjecture C2, they will 
eventually dominate the path integral. Thus, to the extent to which the principle of chiral 
ordering reflects reality (i.e. the extent to which the effective interaction phases affecting 
chiral fermion encapsulate the physics), we suggest that it is the elements ofU L ' D that should 
be examined in order to determine the nature of collective degrees of freedom relevant in 
QCD vacuum. 

(4) Note that the actions reached via the procedure of chiral ordering, while describing 
pure-glue theory, become functions of lattice Dirac kernel D describing the associated chiral 
fermions. This "appearance" of structure associated with fermions in gauge part of the 
action will be a recurring theme in this series of papers and, in fact, it is very natural in the 
context of full QCD. Indeed, we will elaborate on the usefulness of tying the gauge and the 
fermion structure later in this paper and extensively in the in the second paper of this series. 

(5) We emphasize that everything discussed here for the case of exact chiral ordering (Con- 
jecture C2 in particular), also applies for various cases of approximate chiral ordering that 
will be defined in Sec. 5. This is important since it gives the possibility to test the ideas 
proposed here, namely by using the overlap Dirac operator as a starting point. We also 
remind the reader that the set of perfectly chirally ordered configurations U L,D is already 
defined for all chirally symmetric operators D e S F . 

4.5 More on the Rationale for the Principle of Chiral Ordering 

It is sometimes emphasized that the physical content of the gauge field is associated with its 
influence on a charged particle. This is usually expressed by saying that the wave function 
of a particle (charged under the gauge group) acquires a non-integrable (path-dependent) 
matrix phase when traveling over a specific path in space-time [35, 36]. Let us impose this 
condition in the lattice theory for the fermion in the fundamental representation and an 
"infinitesimal" path so that we can avoid the discussion about restricting the particle to a 
particular macroscopic path and/or about the interference of contributions from different 
paths. The elementary displacement of the particle in the lattice theory at finite cutoff is 
not really "infinitesimal" of course. However, we will impose it as if it was in order to mimic 
the behavior in the continuum theory. The quantum dynamics of the particle is described 
by its propagator, and for Dirac particle of mass m we obtain the condition 39 

(D + m)-^(U) = (Df + mtf xU^ Vn,/x,Vf/ U n ,, e SU(N) (74) 

39 It is worth noting that the above motivating discussion in the continuum uses the Hamiltonian language 
for the probing particle. This is actually problematic in the lattice theory with Ginsparg- Wilson fermions 
due to non-ultralocality [34] (see also [37]) and thus the discussion cannot be carried out entirely within the 
lattice context. Thus, Eq. (74) should be viewed as a discretized version of the relation in the continuum 
where the connection of the propagator formalism to Hamiltonian is formally made. 
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where ( + m)~ 1 is a free (translation invariant) propagator (no color indices). The above 
equation is to be compared with Eq. (69). 

Assuming that the fermionic actions satisfying the above condition exist, or using its 
approximate version with arbitrary D G S F , the map M. : U — > U represents a valid 
deformation of the gauge action for non-zero positive m (constant in lattice units). Indeed, 
the only concern is locality which is however expected to be present for finite lattice mass. 
One could thus view Eq. (74) as a condition defining the chiral ordering along the similar lines 
as we did for standard chiral ordering Eq. (69). However, the above line of reasoning would 
not have much meaning if the configuration U depended strongly on the mass m chosen in 
the map. In other words, it would not be natural if the intended definition of physical content 
of U depended strongly on the mass of the probing particle. The consequences associated 
with demanding the independence of map M. on m can be analyzed at large m. Indeed, the 
spectrum of D G S F is bounded, and we can choose the value of m larger than its upper 
bound, i.e. m > 2 for the standard overlap operator. We then have 

x 1 1 / D D 2 D 3 x 
(D + m)- 1 = -(!--+ +...) 

which is a converging Taylor-like expansion in t = 1/m. The requirement that the configu- 
ration U is independent of m then directly leads to the set of conditions 

D n>n+lt (U) = D f x U n ^ , (D 2 ) n , n+ ,(U) = {D% x , ... (75) 

for arbitrary positive powers of D. Thus the chiral ordering condition Eq. (69) represents 
a leading term (in t) which facilitates the independence of the "propagating phase" on 
the mass of the (heavy) probing fermion. This further supports our proposition that the 
effective interaction phase associated with elementary hopping of chiral lattice fermion tends 
to preserve the lattice physics. 

We finally add two relevant remarks. 

(1) From our discussion in this section it is tempting to conclude that the chiral ordering 
prescription given by Eq. (74) is more fundamental than that of Eq. (69), and that the 
principle of chiral ordering formulated in Sec. 4.2 should be taken with respect to this relation. 
This may in fact turn out to be true. However, from Eq. (75) we see that the two are quite 
closely related. Moreover, the perfectly ordered configurations with respect to (74) form a 
subset of those that are perfect under (69). We will thus not loose any physical information by 
proceeding with the original principle of chiral ordering which can be numerically investigated 
with computational power currently available. 

(2) Classical configurations (see footnote 9) belong to the set of perfectly chirally ordered 
configurations in the classical continuum limit, i.e. on an asymptotically large lattice. This 
is true for both versions based on Eqs. (69) and (74). To make it more precise, consider 
a classical field A^(x) on a symmetric torus of physical size L p . Introducing a "classical" 
lattice spacing via L p = La we construct the corresponding configuration U c (A ll ,a) on an 
L 4 lattice via Eq. (9) (see also footnote 10). We then have in case of Eqs. (69) that 

D n , n+ ,(U%A^, a)) - Dl x U%A^, a) = 0(a 2 ) (76) 
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i.e. the linear term in a vanishes. In this classical case it is a simple consequence of gauge 
invariance. We emphasize this fact to reassure ourselves that the principle of chiral ordering 
is compatible with the requirement of classical correspondence. It obviously doesn't imply 
that classical configurations necessarily dominate the path integral for actions at fixed 
large j in the quantum continuum limit. 



5 Approximate Chiral Ordering 

Our discussion of the principle of chiral ordering in the previous section was mostly based 
on the assumption that there exist valid fermionic actions (elements of S F ) for which the 
condition (69) is strictly satisfied (elements of S F ' C ). This allowed us to proceed with our 
reasoning in a straightforward manner. However, we also made it clear that the basic idea is 
meaningful (and can be pursued) without this assumption. In this section we define various 
forms of approximate chiral ordering procedure (map M. D ) which can be realized e.g. with 
the standard overlap operator [12]. It will turn out that our discussion will also shed some 
light on how should one view the possible existence of (non-empty) set S F ' C . 



5.1 Various Forms of Approximate Chiral Ordering 

The proposition that the physical content of the gauge field is largely associated with the 
interaction phase of the probing chiral fermion can be implemented in many different ways. 
In an ideal case, these would all lead to equivalent behavior. However, it is also possible 
that some are superior to others in the sense that the degree to which the procedure tends 
to preserve the physical content of the configuration/theory can vary. We discuss some of 
the possibilities bellow. If not stated explicitly, the maps defined bellow are understood to 
be elements of £F, i.e. they satisfy conditions (i)-(iii) of Sec. 4.1. 

(I) The simplest way to proceed is to consider the formula (69) specifying the chirally ordered 
link U n ,n m the exact case, and use it as a basis for the ordering map also in the inexact 
case. The resulting map M D is then a composition of three maps M D = M.^ oM.^ o 
specified by the following steps. 

(a) Projection to color space, i.e. 



M {1) : U, 



M = 



4 tr * 



(U) 



dI x r)- 1 D 



(U) 



Dl)' 1 D n , n+ ,(U) 



(77) 



where tr s denotes a spinorial trace, 40 I c is an identity in the color space, and I is a configu- 
ration of such identity matrices. 41 The last form is just a shorthand for the first two in this 

40 If O is a matrix in color-spin space, then tr s O is the matrix in the color space with matrix elements 
(tr s 0) a .b = J2 a O aa ,ab, where a denotes the spinorial and a, b the color indices. 

41 Note that if D' happens to be non-invertible (which can occur e.g. in the case of the overlap operator 
for special values of negative mass and the Wilson parameter), the map M (1) is still well defined by Eq. (69). 
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definition. One can easily see that in the case of exact chiral ordering (69) we would have 
U ntli = M and thus M D = M {1) . 

(b) Unitary projection. This is provided by the polar decomposition which offers a natural 
and unique definition of the "matrix phase" for a non-singular square matrix in direct analogy 
to the phase of a non-zero complex number. According to the polar decomposition theorem, 
if matrix M is non-singular then it can be uniquely decomposed as 

M = M u M h such that (M^M" = I c , (M /l ) t = M h (78) 

and M h is positive semidefinite. This associates a "matrix phase" (unitary matrix M u ) with 
matrix M. 42 Thus the map Ai^ is defined by 43 

M {2) : M — > M u = M ; 1 (79) 

where denotes a positive semidefinite branch of the square-root. 44 

(c) SU(N) projection. Unitary matrix M u can be mapped on the element of SU(N) via 
rescaling by the N-th root of its determinant. Selecting the N-th root that leads to minimal 
modification of M u (in arbitrary matrix norm) specifies a unique prescription, namely 

M (3) . M u > jj^ = e ~i<p/N M u det M u = e i<p ^ e ^ (gg) 

The full map Ai D = Ai^ o Ai^ o Ai^ can be expressed in a single explicit formula if so 
desired, and is an element of H. 

(II) Using the same logic as in method (I), we can choose to perform the approximate 
chiral ordering via more "continuum-like" formula than that of Eq. (77). Indeed, consider 
a Clifford decomposition of namely 

16 

Dl = Y,B;r a , a = l,2,...16 (81) 

a=l 

where B® are complex numbers and T a are the elements of an orthogonal Clifford basis 
normalized such that trT a r 6 = A5 a> b- Let us denote by the Clifford component of 
corresponding to basis element 7„, i.e. 

B^ = Bl such that T a = 7l/ (82) 
We then define the projection to color space (map AA^) via 

M (1) : — >■ M = ^tr s [( 7M x I c ) D n , n+ »{U) } = -^tr^D^U) (83) 



IB, 



42 Note that the decomposition also holds for singular matrices but in that case it is not unique. Also, the 
matrix phase M u defined by (78) should be called the left matrix phase since one can also define the right 
phase by switching the order of matrices. Left and right phases are equal for normal matrices M. 

43 Projection to the gauge group based on the polar decomposition (steps (b) and (c) here) was used by 
several authors in different contexts. Reference [38] is the earliest we are aware of. 

44 Wc implicitly assume here that the cases when M is strictly singular will essentially never be encountered 
in practice. One can also take more strict attitude and use the polar decomposition with an additional 
condition to fix M u uniquely also in singular cases. We will not elaborate on this here. 
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where the second form in the above equation is an abbreviation for the first form. One can 
easily check that if D e S F,C (exact chiral ordering) we would have U ntfl = M and thus 
M D = M^. In the approximate case we define M D — o o M^ l \ where 
and are specified in method (I) above. 

(Ill) The non-trivial aspect of the chiral ordering condition (69) resides in the fact that it 
forces a factorization of local spin-color matrix in the respective degrees of freedom. This 
gives meaning to the statement that the fermion acquires a matrix phase when it hops from 
site n to site n + /i. In other words, it ensures the fact that all fermionic degrees of freedom 
at site n get rotated by the gauge field in an identical manner. It is thus sensible to define 
the approximate chiral ordering prescription based on the "degree of factorizability" . For 
example, we can construct the map Ai in the same way as in case I but replace M.^ with 

M {1) : U n „ — >M such that min F(M) — F(M) (84) 

MeC NxN 

where <C NxN is the set of N x N complex matrices and 

F(M) = HA^+^-^xMH (85) 

with ||.|| denoting a matrix norm. One can easily check that for non-zero the solution 
of the above minimization problem in Frobenius norm (i.e. ||^4|| 2 = trA^A for matrix A), is 
unique and given by 

1 



M = - f r tr s 



(Dff D n>n+fl (U) (86) 



Comparing the above expression with Eq. (77) we see that the two projections to color space 
are identical if 

(DffDl = tr K D A tD A V (87) 

i.e. when is a unitary matrix rescaled by a real number. This is true e.g. for the overlap 
Dirac operator but not necessarily in general. Note that if D e S F,C (case of exact chiral 
ordering) then it has to satisfy both Eqs. (77,86) and, consequently, the unitarity condition 
(87) also has to be satisfied. 

Let us also point out that when replacing Dj[ in the minimization problem (84,85) with 
arbitrary spinorial matrix S, then the solution is given by (86) with replaced by S. Con- 
sequently, one can easily check that the color-space projection Ai^ of case II is equivalent 
to minimization of F(M) for S = (fi not summed). 

If one wishes to focus the definition of color projection on the spin-color factorizability 
property, then one can obtain a prescription that is different (but approximately equal in 
generic case) from those in cases I and II by relaxing the condition that the spinorial factor 
is fixed. More precisely, one can associate the projection to the color space with the following 
minimization problem. 

M {1} : U nix — >M such that min F 1 (M,S)=F 1 (M,S) (88) 

where 

F 1 {M,S) ee \\D n , n+ ^(U)-SxM\\ (89) 
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Note that the above problem is invariant under simultaneous rescaling S — > aS, M — > a _1 M, 
and we fix the resulting degeneracy of solutions by requiring that the 7 M -component in 
Clifford expansion of S coincides with that of Dj[, i.e. 

tr^S = AB^ (90) 

At the numerical level the solution can be obtained e.g. iteratively by alternating the solu- 
tions of type (86) with spinorial/color matrix factors fixed. The resulting map A4^ is then 
composed with and of I to obtain the full map M D . 

(IV) In the remaining three cases we list prescriptions that are based on the spin-color 
factorizability property of exact chiral ordering, but avoid steps (b), (c) (i.e. unitarity 
projection and SU(N) projection). In other words, they are based on a direct minimization 
in SU(N) (rather than C NxN ). As a result, they will not be expressed in the explicit form, but 
rather defined only as solutions to the corresponding minimization problems. The uniqueness 
is hard to demonstrate in these cases, and will be assumed. The actual numerical work will 
reveal if this is reasonable or if it will be necessary to impose additional conditions to make 
them unique. We thus define the map M. D for the current case via the solution to the 
following minimization problem 



M D : U n>li — > U n>li such that y min F(V) = F(U n „) (91) 



where 



F(V) = \\D ntn+ ,(U)-DlxV\\ (92) 

The above map represents perhaps the most natural generalization of exact chiral ordering 
prescription (69) to arbitrary D G S F . 

(V) The prescription that formally looks more "continuum-like" and represents the analog 
of case II is given by map (91) with F(V) given by 

F(V) ee 1 1 D n , n+lx {U) - B^lv x V\\ (93) 

(VI) Finally, the analog of case III (Eqs. (88-90)) is specified by the minimization problem 
M D : U n ,„ — > U n ,„ such that min Fi(V,S) — Fi(U n fi ,S) (94) 

VGSl/(Af),5GC 4 >< 4 

where 

F 1 (V,S) ee \\D n , n+ ^U)-SxV\\ (95) 

Note that since V is constrained to be unitary, it is not necessary to impose the analog of 
condition (90) in this case for SU(3) group. 

5.2 Repeated Approximate Chiral Ordering 

Various versions of approximate chiral ordering defined in the previous section represent some 
of the natural choices but certainly do not exhaust all possibilities. In fact, given arbitrary 
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D E S F , we have not explicitly denned the set of acceptable maps M D that we would 
consider approximate chiral orderings. Let us thus put forward the following definition. 

Definition 2. Let D{U) G S F be an arbitrary chiral lattice Dirac operator. We call the map 
Ai D an approximate chiral ordering defined by D if (a) M. D e J '• (b) the elements oftl L ' D 
(perfectly ordered configurations) are its fixed points, namely 

M D {U) = U , U e U L ' D (96) 

(c) it approximately satisfies the chiral ordering relation for generic configurations, i.e. 

e _ lUWtQ - oj> 'AQOII ;<1 v m 



As usual, by "generic" we mean all up to a possible set of measure zero which we assume 
not to be relevant. 

The reason for introducing the approximate chiral ordering is that basically all the fea- 
tures discussed for the exact case are expected to carry over generically to the approximate 
case. In other words, the principle of chiral ordering (see Sec. 4.2) is viewed as a robust 
statement about the nature of the gauge field. In that regard, we wish to explicitly mention 
the following points. 

(1) As emphasized in Sec. 4.3, we assume that there exist approximate chiral orderings for 
which the ordering conjecture CI2 is valid. 

(2) Fixing arbitrary D e S F and a particular approximate chiral ordering .M D , we can con- 
sider the sequences { U { i\ j = 0, 1, 2, . . . } and { S$ D , j = 0, 1, 2, . . . } defined by repeated 
application of M. D as in the exact case. As j increases, it is expected that the configurations 
(and the ensembles representing the theory) will become less affected by randomness with 
lattice physics changing very mildly. Very important expectation is that the evolved config- 
urations will asymptotically approach the elements of U L,D and the ensembles will approach 
a non-trivial lattice theory. To express this, we modify the conjecture C2 of exact chiral 
ordering as follows. 

Conjecture C2a. Let U G U L and D e S F . There exist approximate chiral orderings 
Ai D such that (a) lim^oo exists generically and is an element ofti L ' D ; (b) there exist 
configurations U for which this limit is non-trivial at sufficiently large L. 

Needless to say, performing numerical tests of the principle of chiral ordering (as well 
as conjectures CI and C2a) in the context of overlap Dirac operator will be a crucial step 
towards making the construction proposed here viable. 



5.3 Dual View of the Chiral Ordering Evolution 

Let us recall that the evolution of arbitrary configuration U under chiral ordering (exact 
or approximate) corresponding to D £ S F is characterized by the sequence { j = 
0, 1, 2, . . . }, which in turn induces the evolution of gauge theory S in the set of actions, namely 
{ Sq[) ,j = 0,1,2,...} where the actions are defined via the associated configuration-based 
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deformation. However, it is also sometimes useful to think of chiral ordering as an evolution of 
the underlying chiral fermionic action or equivalently a Dirac operator. Let us thus consider 
a sequence of operators { D^\j = 0, 1, 2, ... } defined via 

D^{U) = D(U&) = D([M D y(U)) j = 1,2,... D<> \U) = D(U) (98) 

Since A4 D is an element of £F, it is easy to check that e S F for arbitrary j. Moreover, 
the free Dirac operators associated with coincide with that of D, i.e. 

D^(l) = D(l) j = l,2,... (99) 

The basic property of built in by construction is that 

= (M D ) j (U) = M DU ~ 1] (U) (100) 

In other words, many iterations of chiral ordering using D can be thought of as a single 
iteration with another element of S F . Assuming the validity of conjecture C2a (C2) we can 
conclude that for generic configurations U there exist limits 

lim U ij) = U lim D {j \U) = D(U) (101) 

satisfying the relations 

DnAU) = D n ^U) = D n ^{U) = Dl x U n „ (102) 

We can then conclude the following, (i) The operator D satisfies the exact chiral ordering 
property (69) regardless of whether the initial operator D satisfies it or not. In fact, it 
is an operator which yields a fixed-point configuration (perfectly ordered configuration; 
element of U L,D ) in one step, (ii) Operators D and D share some of their perfectly ordered 
configurations. 

The above considerations suggest that the validity of conjecture C2a is sufficient to 
conclude that the set S F ' C C S F is non-empty. Indeed, every of chirally ordered 
sequence is an element of S F , and one thus naturally expects that the limit D belongs there 
as well, in which case D G S F ' C . However, we caution that the limiting action could become 
asymptotically non-local in some cases. Nevertheless, even if that happens, the important 
result of these considerations is that the elements of S F can satisfy the chiral ordering 
condition (69) to arbitrary precision. 

6 Partial Chiral Ordering 

In our treatment of the principle of chiral ordering we required that the ordering condition 
(69) be satisfied for all lattice links. In other words, well-defined interaction phases were 
assumed to exist for all possible elementary space-time "hops" of chiral fermion, and for all 
(generic) configurations. However, one could also take a view (and some results of [5] could 
be interpreted that way) that there is a sizable fraction of elementary space-time filaments 
(links) which are not used by the propagating fermion. For such links the interaction phase 
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would not have to be well-defined. Turning the argument around, one could also say that 
the regions of space-time where condition (69) is satisfied to high accuracy are highly sus- 
ceptible to fermion propagation, while the ones where (69) is satisfied poorly form sort of an 
"excluded" space-time. 

To incorporate the above possibility into our reasoning (and to be able to put it eventually 
to the test), requires few changes in the formalism and we will describe them in this section. 
However, our overall line of logic will not change at all. In fact, we will essentially trace the 
same steps as we did in case of full chiral ordering, and will mostly concentrate on what 
needs modification. Note also that we will not use different notation for analogous concepts 
in the two cases which will always be distinguished by the context (either "full" or "partial" 
chiral ordering) in the future. 

(%) To build the concept of chiral ordering we started with the definition of the special class 
of fermionic actions S F ' C , which in the case of full chiral ordering was defined as the subset 
of elements from S F satisfying condition (69). To modify this for the partial case let Q\ be 
the set of all link coordinates (i.e. all (n, /i) of a given lattice) and L C Qi its arbitrary 
subset. We say that D £ S F ' C if D £ S F and if for generic configuration U there exists a 
non-empty L(U) such that 

D n , n+lx {U) = Dl x V (n, fi) £ L U n>li £ SU(N) (103) 

The associated partial chiral ordering map M D is defined via 

Thus, the gauge links for which the corresponding interaction phase is not well defined will 
not get changed under the transformation. One can easily see that MP £ £F thus inducing 
a valid deformation of the gauge action. It is expected that the chiral ordering conjecture 
CI2 applies and the Kolmogorov entropy of S M will not grow relative to S. 

(ii) The analysis of the repeated exact partial chiral ordering has the same ingredients as in 
the exact case. For arbitrary D £ S F we can define the set of perfectly ordered configurations 
U L,D in complete analogy with Definition 1 with Eq. (73) required to be valid only on some 
non-empty set £j(U) of links. Direct analogues of conjectures CI and C2 are expected to 
hold. 

(Hi) The implementation of approximate partial chiral ordering requires a non-trivial modi- 
fication relative to the full case. Indeed, it is not possible just to take e.g. the cases (I)-(VI) 
of Sec. 5.1 and restrict them to a subset of links. The point is that for generic D £ S F 
the equation (103) will not be satisfied exactly for any links, and thus the selection of L is 
part of the constructed map. The set L is configuration-dependent and the rule specifying 
it must be chosen such that the resulting Ai D is local. Below we discuss few examples that 
conform to our requirements. 

(Ip) Consider the map M. D for case (IV) of full approximate chiral ordering with arbitrary 
D £ S F . For given U £ IA L let e nifi (U) be the relative residues of Eq. (97), and let 5 > 0. 
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Then we define the map Mf by 

U — Mf(U) (Mf) n AU) = | "ZJU) ' < * (105) 

We emphasize that .Mf is local and thus the element of 5F. Depending on the value of 5 it 
will chirally order all the links or some subset L. 45 We can obviously base the construction 
of Aif on other full approximate chiral orderings in an analogous manner. 

(Hp) Qualitatively different kind of prescription can be obtained as follows. Consider the 
matrix element D n n+fM (U) and the configuration U(V) obtained from U by replacing U nifl 
with V G SU(N) while keeping all the other link variables the same i.e. 

U mv (V) = \Y> M = {n^) 

, V ; \U m ,v, {m,u)^(n,n) v ; 

Furthermore, let 2J njM (£/) denote the set of all V G SU(N) such that D n:fl+fM (U(V)) leads to 
an exact spin-color factorization of chiral ordering type, i.e. 

W n jU) = {Ve SU(N) : D n<n+ ,(U(V)) =D£x U n ^(V) , U n ^(V) G SU(N) } (107) 

and let Vq be the element of 9Jn jjU closest to U njfl , i.e. such that 

mm\\U n ^-V\\ = \\U n ^-V \\ (108) 

In the last equation we have implicitly assumed that if 2J„ iM is non-empty then V is unique 
generically. We then define 

The above equation defines the approximate chiral ordering Ai D (element of £F) with the 
property that it only modifies links for which the corresponding spin-color factorization is 
achievable in a prescribed manner. It is thus a partial chiral ordering prescription. 

(IIIp) The disadvantage of the partial chiral ordering defined above is that it is not clear 
at this point how to carry it out in practice. In other words, there are no straightforward 
numerical ways to determine Vo, and they will have to be developed. We thus also advocate 
another possibility. Recall that one of the main purposes of approximate chiral ordering is 
to determine the structure of perfectly chirally ordered configurations obtained by repeated 
application of map Ai D . A natural way to proceed in this regard is to enforce "local perfect- 
ness" in the map Ai D . By this we mean the following. Using the same notation as in case 
(Hp) let 2J£ ,u(U) denotes the set of all possible "perfect links" U n ^ associated with fixed 
(n,/i), i.e. 

% P njU) = {U^ G SU(N) : D n , n+ ,(u(U n , p )) = D[ x U n ^} (110) 

45 Note that S is not to be viewed as a free parameter that one can change at will when changing the 
configuration or the lattice spacing. The map Mf has to be fixed when investigating the evolution in the 
set of actions. 
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We now select the element U® ^ of QJ£ ,«(^0 an y) closest to U n ^, i.e. such that 



minll^-^H = \\U n> ,-U^\\ (111) 
and define the map Ai D via 

_ ItW, if 93^ = 1 } 

Note that the important difference between the sets 2J„ iM (£7) and 9J£ „(£/) is that in the 
former case we have twice as many variables to solve for defining condition than in the latter 
case (16 versus 8 in SU(3) case), and thus the set 03^ ,iSP) c %3n,n(U) is naively expected to 
be quite small generically. In fact, one expects at most a finite set of solutions in the perfect 
case. However, the comparison to case (Hp) can only be made by studying these solutions 
on specific backgrounds typical for the starting theory S. Practically relevant feature of the 
perfect case is that one can attempt to find the elements of %S^^{U) in a straightforward 
manner iteratively. 

(iv) The content of conjecture C2a carries over naturally to the partial chiral ordering case 
and the analogous statement is expected to be valid. One can also take a dual view and 
think in terms of evolving Dirac operator as in the full case. 

Let us finally remark that the usefulness of the idea of partial chiral ordering can only 
be decided by studying the above options (full and partial) explicitly. Indeed, it is an 
open question whether QCD dynamics leads to the phenomenon of excluded space-time, i.e. 
whether the natural splitting of space-time into regions of well-defined and ill-defined chiral 
phase actually occurs in typical configurations dominating the QCD path integral. 



7 Chiral Ordering and the Renormalization Group 

As discussed extensively in the previous parts of this paper, the search for the fundamen- 
tal QCD vacuum structure in the path integral formalism naturally leads to exploring the 
behavior of Kolmogorov entropy in the set of lattice QCD actions. Moreover, we proposed 
that a suitable tool to guide us in this search might be the iterated chiral ordering transfor- 
mations defining certain trajectories in this set. Thus, at least superficially, there exist some 
similarities here with the Renormalization Group (RG) concepts and techniques. The aim 
of this section is to suggest, and make it plausible, that the connection might be more than 
superficial. 

In case of RG the important structure in the set of actions is represented by renormalized 
trajectories associated with different RG transformations. The lattice actions on renormal- 
ized trajectory are expected to yield exact predictions (no cutoff effects) for all spectral 
quantities (or any physical quantities if the operators are correspondingly improved) [39]. 
While not clear a priori at all, it is tempting to think that renormalized trajectories will also 
represent a relevant feature from the point of view of Kolmogorov entropy. Such expectation 
is tied to the fact that, in certain situations, it appears valid to intuitively associate the phys- 
ical (continuum) content of the lattice theory with order in typical configurations (clearly 
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defined collective variable C), while to view scaling violations as being caused by an excess of 
"disorder" superimposed on top of it. Indeed, consider the theory S G (a({3)) G S G and follow 
the trajectory of actions as a is being changed towards the continuum limit. As the lattice 
spacing is decreased, the scaling violations become smaller while the Kolmogorov entropy 
is expected to decrease since the associated distribution of configurations is becoming more 
sharply peaked. In this sense the transition from lattice spacing a\ to lattice spacing a 2 < ai 
can be intuitively viewed as a "noise reduction" in typical configurations of S G (a\). Looking 
now in the orthogonal direction in the set of actions (lattice spacing fixed), the actions on 
renormalized trajectories represent the solutions to the problem of maximizing the infor- 
mation content of the lattice theory relative to the continuum. In other words, among the 
actions at the particular fixed cutoff, they carry the most information about the continuum 
limit. As such, the underlying order shaping this physics is expected to be closely related to 
the one associated with the physics of the continuum limit. Thus, in analogy to approaching 
the continuum limit, it is natural to expect that departures from this specific order (moving 
away from renormalized trajectory) are caused by extra randomness present in the dynamics 
defined by the corresponding actions. In other words, one tends to conclude that actions on 
renormalized trajectories correspond to local minima of Kolmogorov entropy at fixed lattice 
spacing. 

While the above considerations appear quite plausible, they are strictly heuristic. At the 
same time, the conclusion on renormalized trajectories is rather non-trivial. 46 In this section 
we extend the concepts needed to introduce the principle of chiral ordering in a manner that 
brings the relation to the renormalization group ideas more to the forefront. We will then 
argue that, in this extended framework, the above conclusion on renormalized trajectories 
suggests itself naturally. In fact, we will propose that the iterated chiral ordering can drive 
gauge actions close to renormalized trajectories of particular RG transformations. 



7.1 RG Transformations Based on Chiral Ordering 

Let D G S F and Ai D be an arbitrary full approximate chiral ordering map. 47 The action 
of M. D on U ntfl is determined by D ntn+lx (U) and it represents the interaction matrix phase 
associated with "elementary hopping" of chiral fermion from n + ji to n. In the same manner, 
we can think of an effective interaction matrix phase associated with fermion moving from 
n + 2/i to n. Such map will be based on -D nj „+2 M but otherwise be identical (in form) to M. D . 
To illustrate this explicitly, let us consider the approximate chiral ordering Ai D of case (IV) 
(see Eqs. (91,92)). The associated map M D ' 2 will be given by 

M D > 2 : — > M%(U) = U n ,„ such that y mm^ F(V) = F(U n ^) (113) 

46 Wc point out that in case of approaching the continuum limit in given lattice theory, one can support 
the conclusion of minimal Kolmogorov entropy with the fact that all correlation lengths increase sufficiently 
close to the limit. This is not the case when approaching the action on renormalized trajectory in which case 
the correlation lengths are merely adjusted to reproduce the continuum masses. This is one of the reasons 
why the conclusion of minimal Kolmogorov entropy is non-trivial in this case. 

47 Note that if it happens that D G S F ' C then all approximate chiral ordering prescriptions define the same 
map, namely the associated exact chiral ordering. 
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where 

F(V) = \\D n ^{U)-D^xV\\ (114) 

In general, we associate with arbitrary approximate chiral ordering map M. D the family of 
maps M D ' S defined by the formal replacement 

M D ^ M D,s under \Dn,n + ,(U) -^D^ n+sfx (U) g = ^ 2> . . . (H 5 ) 

[ H > SP- 

where obviously M. 0,1 = M. D . 

It is now straightforward to define the associated family of RG transformations with 
scaling factor s. For simplicity, we will work formally on the infinite lattice. The discussion 
can be repeated for the finite lattice in a straightforward manner if one lets the lattice sizes 
L to be the appropriate factors of powers of s. The RG transformation is defined by 

U _> ^{U) : U nix _> $D£(y) = M%{U) (116) 

where reverting to the original lattice is already performed. Note that "Jl 0,1 = Ai D . Repeated 
application of 3l D ' s on U induces the sequence of configurations U s ^ = (Oi D,s y(U), j = 
0,1,2,... with U s, ° = U. Starting from the gauge action S = S^' s G S G represented by 

the corresponding ensemble, this repeated map defines the sequence of actions S^' s , j = 
0,1,2,... in S G , namely the RG trajectory. In fact, because of the freedom in choosing s, 

we have defined a family of RG trajectories with different scaling factors s = 2, 3, Chiral 

ordering procedure (s = 1) represents the limiting case in this family and we will refer to 
the corresponding trajectory as the chiral ordering trajectory. 

To close this section, we need to discuss a particular aspect of the above transformations 
that one needs to keep in mind. In a generic case of RG transformation with scaling factor s, 
one argues on intuitive physical grounds that the long-distance behavior in the induced theory 
of blocked variables is exactly the same as in the original theory. This is expressed in the fact 
that the correlation lengths (before reverting to the original lattice) remain exactly the same. 
As a result we find (after reverting to the original lattice) that — > ^/s in the transformed 
theory for all correlation lengths £j. Consequently, the lattice spacing (being based on one 
of the correlation lengths) gets magnified as a — > sa. The conclusion that interactions on 
renormalized trajectory specify perfect actions (no cutoff effects) relies heavily on the fact 
that the above scaling of correlation lengths is exact. This is difficult to prove rigorously for 
generic non-linear RG transformations. Moreover, even in case of linear transformations (e.g. 
for scalar fields) the argument assumes that the transformations are strictly ultralocal, i.e. 
that the new variable does not depend at all on old variables outside the considered block. At 
the same time, it is not difficult to see that if the dependence on variables outside the block 
is non-ultralocal but local, i.e. that the influence of variables outside the block is bounded 
by the exponential decay with characteristic lattice length a, then the correlation lengths 
of blocked variables can in principle change at most by a. If such change occurs, it will 
certainly make little difference close to the critical surface (g — a — 0) when the correlation 
lengths are very large. However, in principle it can make a difference when transforming a 
system with lattice correlation lengths of order a. 
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The RG transformations proposed here are based on Ginsparg-Wilson fermionic ker- 
nels that are necessarily non-ultralocal in fermionic variables [34] and expected to be non- 
ultralocal also in gauge variables which is of relevance here. Consequently, we will not claim 
that actions on the corresponding renormalized trajectories are perfect at arbitrary lattice 
spacing (correlation length). Rather, we will implicitly understand that they are close to 
perfect down to correlation lengths comparable to a. 48 Since there are Ginsparg-Wilson 
actions with a ~ 1 (see e.g. Ref. [40]), this is not an essential restriction. Moreover, as 
we emphasized all along, the nature of the RG transformation based on chiral ordering is 
expected to be such that the physics of the gauge field remains largely preserved even at 
small lattice distances. We thus expect that the correlation lengths will only change by 
the amount smaller than a. 49 The main purpose of the above remarks is to convey two 
messages. (1) Even though the RG transformations proposed here are non-ultralocal, we 
expect that "perfectness" on the renormalized trajectory is still a valid concept. (2) Because 
of non-ultralocality we admit the possibility that the correlation lengths of blocked variables 
can in principle change by very small amounts. 

7.2 Renormalized Trajectory and the Line of Perfect Chiral Order 

While the chiral ordering transformation M. D and the RG transformations of Eq. (116) have 
been put on the same footing via the family of maps 3l D ' s , s = 1, 2, . . ., we would conven- 
tionally tend to think of them in rather different terms. Starting with RG transformations 
3l D ' s , s = 2, 3, . . ., we accept the standard picture of the RG flow for Yang-Mills theory in 
the set of actions. In other words, for given s we assume that there exists a unique fixed point 
(FP) of "3l D,s on the critical surface (g — a — 0), and that there is a single weakly relevant 
coupling associated with this FP, namely the gauge coupling. Starting with arbitrary action 
S(a(g)) e S G close to critical surface (g very small), the RG trajectory will first run towards 
the position of the FP and then away towards larger couplings (larger lattice spacings). The 
limiting procedure where the bare gauge coupling of the starting action S(a) approaches zero 
(critical surface) defines a trajectory running from FP towards larger couplings, namely the 
renormalized trajectory corresponding to 3l D ' s . The generic flow of actions (starting from S 
not necessarily very close to critical surface) in relation to renormalized trajectory is shown 
schematically in Fig. 1 (left). 

Turning now to the chiral ordering map Ai D , one could naively view it as nothing more 
than a very sophisticated way of smoothing the gauge field. Generic procedures of this 
type are sometimes distrusted if used inappropriately. However, as we argued extensively, 
our suggestion here is that chiral ordering is considerably more fine-tuned than generic 

48 Wc emphasize that we do not claim that non-ultralocality necessarily forces the small change of correla- 
tion length in the transformed system. We only say that this is possible with non-ultralocal transformations. 
Also, if one wants to be more precise, our underlying assumption here is that performing the blocking oper- 
ation arbitrarily many times (without reverting to original lattice) leads asymptotically only to a finite total 
change of correlation length. Then the concept of perfect action is meaningful only up to correlation length 
comparable to this total change which is expected to be of order a. 

49 It is worth mentioning here that standard RG transformations in momentum space correspond to non- 
ultralocal (and sometimes even non-local) transformations in position space. Nevertheless, they are widely 
used albeit for analyzing different questions than those relevant here. 
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a'g) a(g) 

Figure 1: The expected flow under transformations R D,S for s > 1 (left) and s = 1 (right). 
Ki collectively denotes the set of couplings other than the gauge coupling. The line of perfect 
chiral order plays the role analogous to that of the renormalized trajectory. 

smoothing 50 since it tends to preserve the defining feature of the gauge field, namely its 
influence on the charged particle. This point of view is expressed by the principle of chiral 
ordering and, in more definite terms, by conjectures CI, C2 and C2a. These imply that 
the evolution in the set of actions under M. D could actually be very constrained, and the 
tendency to preserve physics leads to the fact that the actions Sfy 13 = S^' 1 ,j = 0, 1, 2, . . . 

converge to the limiting action S(S, M. D ) (fixed point) under repeated application of M. D . 
Fixing the starting lattice theory S(a) parametrized by its lattice spacing, we obtain a one- 
parameter line of fixed points S(a,Ai D ). We will refer to this line as the line of perfect 
chiral order. Note that the lattice spacing of S(a, Ai D ) is not expected to be exactly a, 
but rather a(a) since correlation lengths can change by small amounts upon the repeated 
chiral ordering. However, because of tendency to preserve physics, we expect that a ~ a 
with increased accuracy as a — > 0. We should also mention that by writing S(a,A4 D ) we 
implicitly assumed that the line of perfect chiral order only depends on transformation M. D 
and not on the starting theory chosen (analogue of the uniqueness of the fixed point in the 
RG case). We thus propose a flow pattern for Ai D = Ol 0,1 as shown in Fig. 1 (right). 

Given two arbitrary and unrelated RG transformations, it is not expected that the posi- 
tions of their corresponding fixed points and RG trajectories are highly correlated. However, 
it is clear that various transformations in the family 0l D,s are strongly related by the fact 
that they are all defined from a single and highly "fine-tuned" object - the chirally sym- 
metric lattice Dirac operator. Moreover, all these transformations are based upon the same 
underlying principle: they express the effective matrix phase associated with the propagating 
chiral fermion. Thus, for example, one would expect that the evolution of the action under 
a single s = 4 transformation CR D ' 4 is quite similar to repeating two s = 2 transformations 
5l D ' 2 . Consequently, one would also expect that the corresponding fixed points and RG tra- 
jectories for different scaling factors are closely positioned in the set of actions. At the same 
time, as made clear by the discussion in this section, RG trajectories for s = 2,3,... are 
the direct analogs of the line of perfect chiral order in the s = 1 case. Consequently, it is 

50 By "generic" smoothing we mean any map that tends to average the gauge fields locally. Examples of 
this could be the APE smearing [41] or the "stout link" operation [42]. 
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natural to conclude that the line of perfect chiral order runs in the vicinity of RG trajectories 
corresponding to transformations 3i D ' s . 

Accepting the above scenario entails two noteworthy implications: 

(i) Repeated chiral ordering brings lattice gauge theory close to particular RG trajectories 
and thus, due to the conjecture CI2, the regions around such RT-s are special in that their 
Kolmogorov complexity is very low. Thus, the framework of chiral ordering is compatible 
with the intuitive reasoning put forward at the beginning of this section. 

(ii) To the extent that the arguments put forward here will turn out to be valid, one can also 
expect that chiral ordering might be capable of not only "preserving the physical content" 
as stipulated by Goal 2 , but also to improve the scaling behavior. Indeed, since theories in 
the immediate vicinity of RT are expected to have small cutoff effects, it can happen that 
the various correlation lengths will readjust themselves in order to mimic more closely the 
continuum behavior. 

Finally, let us mention that we carried out this discussion in the context of em full chiral 
ordering. Adaptation to the case of partial chiral ordering is possible but depends to a large 
extent on which particular ordering is used. We thus postpone the detailed discussion for 
this case to dedicated publications. 

8 Effective Lattice QCD I. 

The final topic we will deal with in this first paper of the series is the construction of the 
framework for studying the effective QCD vacuum structure and its scale dependence (see 
Sec. 1.5 and points A3 and A4). More precisely, here we will only discuss a specific issue (in 
fact a particular version of it) relevant for this problem. The topic will be discussed further 
in the second paper, and in its full generality in the third paper of this series. 

We will examine the problem of defining the theory S\ F associated with S G S G , for 
arbitrary fermionic response scale Ap. The particular approach proposed here is based 
on the approximate chiral ordering, and will not allow us to do it for arbitrary S G S . 
However, with every S G S G we will associate another valid lattice action(s) with the same 
physics content (in the spirit of the principle of chiral ordering), for which it can be done. 
Viewed alternatively, the construction applies for arbitrary S G S G but with restricted set of 
observables, i.e. we will not be able to define the mean values of all lattice operators 0(U) 
in S\ F , but only a subset of it. However, this subset is sufficiently rich so that this is in fact 
a very mild restriction. 

Let S G S G , D G S F and let A4 D G 5F be an arbitrary approximate chiral ordering based 
on D. As discussed in Sees. 4.1,4.2, the map Ai D can be used to associate with S another 
gauge theory S via the transformation of the ensemble 

S s = {... XP-\ U\ U i+ \ ...} — {. . . M ^ 1 - 1 ), M D (U l ), M D (U l+1 ), ...}= E sM n 

(H7) 

To define the action of effective lattice QCD at fermionic response scale via its 

ensemble, we first define the chiral ordering map Ai D,Xp at the corresponding scale in lattice 
units (Xf — Apa). To do that, we will proceed in the same way as in definition of the effective 
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topological charge density in Refs. [4, 6], i.e. we will perform the eigenmode expansion of 
the corresponding matrix element of D. Indeed, the map M D depends on the collection of 
spin-color matrices D ni „ +jU and we can write 

D n , n +n = J2^i X j (i> J n+^ Dt + ,= E ^A j (^ + ,) t (118) 

j'-\\j\<\F 

where ^ is an eigenmode of D with eigenvalue \j and, in what follows, we will always 
assume that the labeling of eigenmodes respects the ordering of eigenvalues by magnitude, 
i.e. if j < k then |Aj| < |Afc|. 51 The map A4 d ' Xf is then defined via a replacement 

M^^Mffi under ^ -> ^ D' — > (119) 

Note that in the above equations A^ G [0, |A| max ], where |A| max is the maximal possible 
magnitude of an eigenvalue. The effective lattice QCD Sjf F associated with S M is then 
specified by its ensemble 

£ s = {...U*-\U\U* + \...} — 

— > {...M D ~ aKF (U l - l ),M DraKF (U l ),M D ~ akF (U l+l ),...} = S sM d (120) 

with the new probability distribution (action Sj^ (U)) given by equations (62, 63). Note 
that a = a(a) is the lattice spacing of S M ° . According to the principle of chiral ordering we 
expect that a(a) ~ a. 

Several important points regarding the definition and meaning of effective lattice QCD now 
need to be discussed. 

(i) It needs to be stressed that, at this point, we have not yet given the precise definition of 
D^ Xp (i.e. of the effective free matrix element) in Eq. (119). At the same time, the definition 
of Mn,'u F requires both D^ n+lx and D^ Xp for all prescriptions of approximate chiral ordering 
we discussed except for the case (VI) . We will rectify this in the following subsection devoted 
to this point. 

(ii) We wish to reiterate that the map M D ' Xp associates with arbitrary U G U L an effective 
chirally ordered configuration U Xp G U L since 

lim U Xf = lim M D ' Xp (U) = M D {U) = U (121) 

H"^|max \f~~ H A| m ax 

i.e. U Xf is an expansion of U (not U). This is, of course, why Sj^ F given by (120) defines an 
effective description of S M ° rather than that of S. However, at the same time, following the 
reasoning in Sec. 4.1, we can view the map Ai D as the means of associating with arbitrary 
local operator O a (U), that we wish to evaluate in S, a new local operator (U) = 
O a (Ai D (U)) with the same classical limit. The map Ai D,Xp then defines an effective operator 
qM d Af ( t0 be va i ua ted in S) via 

0^°' Af (U) = O a (M D ' ahF {U)) (122) 



51 Note that the expansion in the form (118) applies when D is a normal matrix, e.g. for the overlap Dirac 
operator. If D is not normal then we can still use the expansion in the bi -orthogonal left right eigensystcm. 
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Thus, from the point of view of action S, we can define the effective description at fermionic 
response scale for all operators that can be obtained as (U) for some O a . 

(Hi) It is tempting to think that for the definition of D x f n+fM (and thus of M d,Xf ) it would be 
advantageous to utilize the trick used in [4, 6] to define the effective topological density i.e. 
to eigenmode-expand (D — |A| max I) n)n+/i = D ntn+fl rather than D ntJl+fJj . itself. This would 
seemingly give the low-lying modes the largest weight. However, contrary to the case of 
topological density, it turns out that 75-Hermiticity forces the two expansions be identical. 

(iv) We emphasize that the purpose and the meaning of effective lattice QCD is quite different 
from that of the common notion of effective field theory. In the latter, one constructs new 
degrees of freedom and couples them via effective interaction designed so that it takes into 
account the influence of high-frequency fluctuations that have been eliminated from the 
theory. However, in case of effective lattice QCD the nature of field variables is preserved 
(they are still the same gauge fields), and the interaction is changed only to filter out the 
high-frequency fluctuations in a meaningful manner. This effective interaction is non-local 
with the corresponding range being controlled by the scale Ap. The continuum limit of 
the theory S ,Ai? (a(a)) is achieved via bringing the driving theory S(a) to the continuum 
limit (a — > 0) while keeping Ap fixed. The values of measurable quantities in the continuum 
theory S^' D will differ from those of the full continuum theory S® CD , and this difference 
must then be directly ascribed to the influence of high-frequency fluctuations in the gauge 
field. In the language of the typical vacuum configurations, the space-time vacuum structure 
acquires a scale-dependence via Ap, and the relation of this structure to the aforementioned 
differences will teach us how is vacuum structure at different scales related to various physical 
phenomena. In other words, the purpose of effective lattice QCD is to provide us with the 
scale-dependent picture of QCD vacuum. 

(v) While we have carried out the discussion of effective lattice QCD in the context of pure- 
glue QCD, it can be clearly repeated without any change for full QCD. In fact, for given 
full lattice QCD defined by S G G S G and S F G S F , there is an obvious natural choice of 
Dirac operator on which we base our chiral ordering transformation A4 D , namely the D that 
specifies S F . 

(vi) At the technical level, our way of defining effective QCD is similar to the recently 
discussed procedure of Laplacian filtering [21], wherein the authors define a "filtered link" via 
performing the eigenmode expansion of relevant matrix element of the covariant Laplacian. 

8.1 Expansion of the Free Matrix Element 

In this section we come back to the issue of specifying D^ Xp , which is necessary for definition 
of Ai D ' XF with approximate chiral orderings we discussed (except case (VI)). Our goal is 
to preserve the physical meaning of chiral ordering procedure in the effective case. In other 
words, we want that the map A4 D,Xp can still be interpreted as extracting the matrix phase 
acquired by fermion when hopping from n + fi to n relative to the free case. The issue is what 
exactly do we mean by "free case" in effective theory. There are several ways to proceed and 
we will describe them in turn. 
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(a) The most straightforward possibility is to expand the free operator (no color indices 
in this case) 52 in its own eigenmodes x J . The relevant matrix element is 

= D f >X = ]T xicvj ( X i + ,y Df X 3 = "jX j (123) 

j:\u>j\<X F 

While this is a sensible definition satisfying our requirement, it might not be suitable if we 
are interested in (or restricted to) studying the low-energy behavior in finite physical volume 
V p = (La) 4 (fixed as continuum limit is approached). Indeed, with antiperiodic boundary 
conditions in time direction, there will be a non-zero lower bound of the free spectrum 
Mmin(£) such that \uj\ m in{L) / a oc l/(La) converges to a finite physical value A™ m as a — > 0, 
when V p is held fixed. Thus, in this case we will not be able to define the map J\A D ' aAF 
for A F < A™ in . In other words, the effective lattice theory SjQ will not be defined for 
A F < Af n . 53 

(P) The possible issue with definition (a) is that for A F close to A™" 1 there will be an 
increasingly larger mismatch between the number of eigenmodes contributing to D^ Xp and 
D^ n+fX as the continuum limit is approached. This comes about as a result of the fact that we 
force the scale (lattice spacing a) determined in interacting theory also on the free dynamics. 
This is not natural. Indeed, one should rather match the free theory to the interacting one 
by the "number of degrees of freedom", which are measured in this case by the number of 
eigenmodes included in the expansions. In other words, it is meaningful to ask what is an 
effective matrix phase acquired by interacting fermion described by A" eigenmodes relative 
to the dynamics described by A" free eigenmodes. To implement this properly, let us recall 
that, due to 75-Hermiticity, the eigenmodes of D with complex eigenvalues are naturally 
paired, i.e. for each such eigenmode i]P of D there is a mode ip k = with complex- 
conjugated eigenvalue (X k = A*). The unpaired eigenmodes have real eigenvalues and reflect 
the topology of the underlying gauge field. The eigenmodes of are all paired via the 75 
operation. For given fixed U, let N be the number of zeromodes of D(U), and let N p (a) is 
the number of paired modes such that \Xj\ < a > 0. By definition, N p (a) is an even number. 
We then have 

Df* = Dl;X = Y: Xi <»j (Xl + »y (124) 

j<N p (X F ) 

Two points need to be emphasized with regard to the above definition. 

(i) Since there can be degeneracies in the free spectrum, we implicitly assume that the 
ordering of the eigenmodes within the degenerate subspace is fixed by some arbitrary rule, 
and that the 75-conjugated modes follow one another when the sequence {(ouj, x*)->3 — 
1,2,.. .} is formed. 

(ii) Note that there is still seemingly a small mismatch between the number of modes con- 
tributing to D^ n+lx (equal to A" + N p (X F )) and the number of modes contributing to D^ Xp 
(equal to N p (X F )). However, one can easily convince himself that the real modes in fact do 
not contribute to D^ n+fi - 



52 Note that since is translation invariant we frequently used (single index) as a shorthand for 
D^ n n+m , Vn. The meaning of single and double-indexing is thus clear. 

53 Note that this is not necessarily a bad thing since l/(La) is a natural infrared scale for this system. 
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Using the above construction, the map M. D ' Xp (and thus action Sf^) is defined by 
equations (118,119) and (124). To see more explicitly how is this done for some basic forms 
of approximate chiral ordering, consider the case (I). We have M D = M {3) o M {2) o M {1) , 
where only the map Ai^ (projection to color space) will need modification when transiting 
from M D to M d ' Xf . In particular, in place of Eq. (77) we will have 

(U n ,,) a , b ^ M Xp = ^[(Dfxir 1 ^)]^ 

= \ E \ (li»)t {D^Y 1 W n ) a (125) 

where (ipi) a is the 4-component object (color index fixed to a) and Dj l ,Xp is given by (124). 
For approximate chiral ordering of case (II) we will have in a similar manner 

(u n ,,u ^ m x i = ^r F ^[i,D x ^{u)] ab 

= ttL E \ (^)! 1, Ma (126) 

W i:|A 3 |<A F 

where B Xf is the 7^ component in the Clifford decomposition of D^ Xp . We proceed in a 
straightforward manner also in the remaining cases of the approximate chiral ordering. 

(7) The final possibility that we wish to discuss treats D nyn+ll (U) and £) n>ri+M (I) = x F 
in a symmetric manner with respect to the eigenmode expansions performed. In particular, 
we will expand both operators in interacting eigenmodes ipi . To do this, let denotes the 
eigenmode of D(I) = x I c with eigenvalue Uj, i.e. has both the spin and color indices 
and is related to eigenmodes of in an obvious manner. We can write D(I) as 

D(I) = £ X j "j{x j ) ] = E c kuk 2 ^ kl (^ 2 ) t c klM = (^)t£)(I)^ (127) 

j k\,k2 

which can be truncated to define 

D^ n+ ,{1) = E c klM ^ (^ 2 )U (128) 

k\,ki 
|Afcj |,|Afc 2 |<A_F 

The aspect that needs to be discussed with regard to the above definition is that, contrary to 
cases (a), (/3) where spin and color variables are a priori separated, the matrix D Xp n+ Jl) is 
not expected to be exactly expressible as a direct product in the corresponding subspaces. If 
this was possible and we could write it as Df l ' Xp x I c , then we would use this correspondence 
as a definition of D^ Xp in this case. However, as it stands, we need to take an additional step. 
While there are other ways to proceed, 54 we wish to explicitly discuss the approach where, 
instead of defining D^ Xp explicitly and proceeding as in cases (a), (/3), we simply replace 
the factors D^ Xp x I c in various definitions of approximate chiral orderings with D Xp n+fM (I). 

54 0ne immediate possibility would be to define D^ Xp by minimizing the norm 11-0^+^(1)— D^ Xp x 
with respect to D^ Xp . 
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To see this explicitly, let us consider some interesting cases of approximate chiral ordering. 
In case (I) we can replace the projection to color space (see Eq. (77)) directly with 



U n>li — M Xp = itr s 



D X ^ n+ »(I))- 1 D^ n+ll (U)\ (129) 



thus defining the map M. D,Xp by composing it with M.^ and M^. For case (IV) we define 
the map M. D,Xp via the following minimization problem. 

M D + P : U n ,» — > C7 n A - such that y mrn^ F(V) = F^-) (130) 

where 

F(V) = \\D^{U) - £>^(I) (I s x V) || (131) 

where I s is the identity in spinor space. 

Finally, we emphasize that practical utility of cases (a) — (7) is to be decided numerically 
by studying which expansion leads to the most accurately satisfied chiral ordering condition. 



9 Summary 

The purpose of this series of articles is to discuss in detail the set of ideas that, we think, 
can form a basis of a consistent framework for systematic study of QCD vacuum structure 
in the path integral formalism. While the history of using the lattice definition of QCD 
to pursue these issues is rather long, it has not matured into a well-defined subject with 
accepted scope, goals, and a standard collection of tools. Thus, while in the areas such 
as hadron spectroscopy one walks on the firm ground of well-defined quantitative notions, 
the realm of QCD vacuum structure is typically associated with constant stumbling over (at 
best) intuitive concepts such as "typical configuration", "object in the vacuum", "space-time 
structure", etc.. More importantly, the ultimate goal of the effort (the analog of obtaining 
the mass) is the part of the research question here. In other words, solving the problem of 
QCD vacuum means different things to different researchers, even though (vaguely) stated 
goals such as "understanding confinement" and "understanding chiral symmetry breaking" 
are of universal interest. 

Nevertheless, the common aspect of various efforts can be summed up rather succinctly 
as follows. We are trying to replace the QCD path integral ensemble S QCD (defined via 
some regularized limiting procedure, e.g. using lattice gauge theory) with the ensemble £ STR 
which is defined in terms of some collective variables C that arise naturally as a result of 
strong dynamics. The analysis of strong interactions in terms of C is expected to be very 
transparent, and the space-time nature of C will represent the space-time structure of the 
QCD vacuum in the path integral formalism. Stating the problem in this way, there are 
two kinds of issues that need to be addressed for there to be a chance that a systematic 
framework emerges, (i) First, it needs to be clarified what is the relation of £® CD (full 
theory of strong interactions) to £ STR that we seek to find (see PI), (it) Secondly, once 
the underlying goal of the effort is specified, it is necessary to have a set of guidelines/tools 
that will make the inquiries into specific questions possible and streamlined (see P2). Our 
suggestions with regard to (i) and (ii) can be schematically described as follows. 
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(i) Given the existence of the fundamental topological structure observable in regularized 

QCD ensembles [5], we propose that the ensemble £ STR [C] can be regarded as fully 
equivalent to £® CD for all questions of physical interest. In other words, we can basi- 
cally view the transition from S QCD [A] to S STR [C] as a complicated change of variables 
in the QCD path integral. Analogously to the case of topological vacuum, we refer 
to the structure entailed by variables C as a fundamental QCD vacuum structure. 55 
Constructing £ STR will not provide us with direct understanding of various phenomena 
in strong interactions because of scale dependence inherent in field-theory description. 
To achieve such understanding, we propose that it is also necessary to study the ef- 
fective vacuum structure which captures the dependence of vacuum properties on the 
scale of fluctuating elementary fields [4, 6]. The effective structure is specified by the 
ensemble £j[ TR = £ STR [C(A),A] describing the effective QCD at scale A. 56 Thus, in 
summary, the ultimate goal in our approach is to construct the ensembles £ STR and 
£[ TR for all A. 

(ii) The above goal can be systematically pursued using the lattice definition of QCD and, in 

particular, by studying the space-time structure in numerically-generated ensembles 
of finite systems in a Bottom-Up manner (see section 1.4). This means that we will 
attempt to arrive at the definition of £ STR and £^ TR inductively, using the observed 
structure in the regularized ensembles as the only input (see A1-A4). In other words, 
the criterion for selecting the appropriate collective variables C (or C(A)) is identified 
in this approach with actual observability of space-time attributes associated with C 
in the regularized path-integral ensembles. 

Having the basic overall scheme described above in mind, the content of this series of 
papers can be described as an attempt to elevate it to the state where its various elements are 
properly defined to the largest extent, and where it is possible to practically pursue its goals 
using the currently available methods of lattice gauge theory. The present work addresses 
three topics related to these issues. 

(I) When searching for fundamental vacuum structure in a Bottom-Up manner, we can 
take advantage of a freedom that is inherent in the process of regularization. In particular, 
the continuum theory can be defined in many ways using different lattice regularizations 
(theories in the set of valid lattice actions). Since the underlying goal is to identify the 
space-time patterns in configurations dominating the regularized path integral, we naturally 
prefer lattice theories for which these configurations exhibit a high degree of space-time or- 
der. However, such reasoning assumes that there exists a well-defined notion for "degree of 
space-time order" in an arbitrary configuration. We argued that Kolmogorov entropy (Kol- 
mogorov complexity) of binary strings describing coarse-grained configurations provides the 
appropriate quantitative measure. The ensemble average of Kolmogorov entropy associated 

55 Note that discovering the nature of variables C does not necessarily imply the analytic solution of QCD. 
However, if QCD can be solved analytically, then specifying C would clearly mean an important step in that 
direction. 

56 For more precise meaning of "effective" in this case see comment (iv) of Sec. 8. More detailed discussion 
will be given in the third paper of this series. 
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with any theory then defines the ranking of regularizations (at a particular fixed cutoff) by 
the degree of space-time order the corresponding lattice interaction generates. We wish to 
make a few comments regarding these issues. 

(a) The Kolmogorov entropy measure is highly universal, capturing the space-time order 
in essentially any form. There is a price that we have to pay for this generality. Indeed, 
if one wishes to consider Kolmogorov entropy as a practical tool (rather than a conceptual 
construct), then one needs to emphasize that the entropy E K (U,k), while well defined, is 
not computable in general. Also, the use of pseudo-random (deterministic) numbers in 
actual simulations will most likely lead to a systematic error in determining the average 
value E K [S, k] for theory S. Thus, if pseudo-random numbers are used for generating actual 
ensembles, then we implicitly assume that our conclusions are not affected by this error. 

(/3) We included a rather detailed discussion of the information theory aspects in Sec. 3. 
This is not just to introduce the Kolmogorov entropy, but also because we believe that the 
language and methods of information theory are both appropriate and fruitful for this area of 
research. Indeed, there are many relevant questions regarding the vacuum structure that can 
be both properly formulated and potentially solved by using the framework of information 
theory. For instance, regarding the issue at hand, one can attempt to formulate less general 
but computable measures (based on the quantity of information) that would mimic the role 
of Kolmogorov entropy. 

(7) It is also worth mentioning that Kolmogorov entropy and the information viewpoint are 
quite useful for providing intuitive understanding in situations that are otherwise difficult 
to grasp. For example, Kolmogorov entropy of the Wilson gauge action at currently studied 
couplings is probably rather high since no observable structure could be identified in its 
typical configurations when standard ultralocal operators for composite fields were used. 
From the information theory viewpoint it is clear that in order to reach the theory with 
significantly lower Kolmogorov entropy, it is necessary that it is related to Wilson theory 
by a highly complex map. In other words, this theory must be defined by an action that is 
difficult to work with computationally. Thus, the actions possessing low Kolmogorov entropy 
are expected to be very complicated. It is for related reasons that we are able to observe 
a well-defined structure in the configurations of overlap-based topological density (which is 
difficult to compute) but not in the ultralocal definitions which are comparably "simple" . 

(II) In order to identify the space-time nature of collective variables C in regularized config- 
urations most easily, one would prefer to work with lattice theories that maximize the degree 
of space-time order. However, it is not sufficient (even as a conceptual goal) to simply mini- 
mize the Kolmogorov entropy in the set of actions (see Goal 1 and Goal 2 of Sec. 4). Indeed, 
one should rather perform a minimization constrained by the requirement that the physics 
content of the theory remains largely preserved in the process. As a tool for achieving this 
in practice, we propose the configuration-based deformation of the action defined via chiral 
ordering transformation. Chiral ordering transformation associates with an arbitrary gauge 
configuration U a new configuration U such that U n>II represents an effective SU(N) matrix 
phase associated with hopping of chiral fermion from n + n to n. Since the transformation 
is constructed to represent the local physical meaning of the gauge field (local influence on a 
charged particle), it is expected that both short and long distance (in lattice units) proper- 
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ties of the theory will change very little upon the deformation (principle of chiral ordering). 
The repeated application of chiral ordering transformation on a given ensemble induces the 
evolution in the set of actions. We propose that there exist nontrivial configurations on a 
finite lattice that are stable under the chiral ordering transformation (perfectly chirally or- 
dered configurations), and that the evolution in the set of actions can lead to theories where 
such configurations play a dominant role. We wish to mention the following points related 
to these issues. 

(a) Two basic forms of chiral ordering transformations were discussed. The full ordering 
where all links of the lattice are forced to acquire their effective value dictated by chiral 
fermion, and the partial ordering where only a subset of all links is changed. The second 
option was introduced in order to incorporate the possibility that the QCD path integral 
is dominated by configurations containing regions that are not accessible to fermion prop- 
agation. Indeed, if that was the case then the effective interaction phase would not have 
to be well defined for links over which the propagation is prohibited. Detailed numerical 
studies will be necessary to determine which option is closer to that being realized by QCD 
dynamics. 

(f3) Connections to RG ideas suggest themselves rather naturally in our framework, and they 
were discussed in Sec. 7. We defined the family of RG transformations with scale factor s via 
chiral orderings extended over straight paths containing s links. The original chiral ordering 
can then be viewed as a limiting case of these transformations with s — 1. We proposed that 
actions on renormalized trajectories corresponding to generic RG transformations represent 
important features in the landscape of Kolmogorov entropy over the set of actions. The 
particular suggestion is that there are deep local minima of Kolmogorov entropy in the 
vicinity of renormalized trajectories. Assuming that chiral ordering evolution indeed leads 
to fixed points belonging to the set of actions, we argued that this conclusion can be made 
plausible at least for renormalized trajectories corresponding to RG transformations based 
on chiral ordering. 

(Ill) In the Bottom-Up approach, the construction of £f TR (effective structure) is to be 
guided by the space-time behavior in effective QCD (S® CD ), which can be defined by its 
ensemble £% CD '. The ensemble £% CD is in turn to be defined as a continuum limit of lattice- 
regularized ensembles £^ QCD with A fixed. However, such regularizations have not been 
defined yet. Here we suggested a particular form of the first step toward the general definition 
of this kind. Guided by the fruitfulness of bringing in the chiral fermion structure into 
definition of lattice topological field, we proposed lattice theories defined at a given fermionic 
response scale A F . This definition is based on the chiral ordering transformation discussed 
here and, in this regard, we wish to briefly mention the following two points. 

(a) It needs to be emphasized that in our framework the search for the fundamental and 
effective structure is completely unified and streamlined under the umbrella of the principle 
of chiral ordering. Indeed, in a nutshell we are contemplating the following chain of logic. 
Consider an arbitrary lattice action S G S G . If the goal is to examine its ensemble to 
obtain an information on the fundamental structure, we suggest that it is beneficial to 
make a transition to the action S M representing a configuration-based deformation of S 
constructed via chiral ordering map A4 D , since this is expected to lower the Kolmogorov 
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entropy. Next, we propose to consider the effective theories S-jQ D of S M ° , whose definition 
is made possible by the fact that theory became a function of Dirac kernel D. This will 
provide information on the effective structure in a manner fully consistent with fundamental 
structure since liniA F _>A™ ax SjQ = S M . 

(/3) As is quite clear from our discussion, the framework proposed here is particularly natural 
in the context of full QCD. Indeed, in this case the choice of D e S F in chiral ordering 
transformations is fixed by the fact that D is part of the definition of the theory itself. 
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